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CHAPTER ei 


tN EROpUCTLION 


Purpose 


It would seem that the primary duty, and continuing 
responsibility of business managers (whether civilian or 
military), is to strive for the best possible economic 
results from the meeounees currently employed by, or avail- 
able to them. 

Linear programming has rapidly become one of the very 
important new techniques available for the management of 
resources of all kinds. The employment of digital computa- 
tion equipment and appropriate digital computer codes 
enhances the scope of linear programming generally, and 
allows its efficacy to be used and explored more effectively. 

The purpose of this thesis is two-fold; first, to 
demonstrate the application of linear programming to problems 
of the Navy Supply Corps and second, to present a new and 
original pragmatic approximation solution to a classical 
unsolved problem in an area related to linear programming. 

In developing these topics for the thesis, it was 
Geemed wise to include a rather full explanation of linear 
programming methodology. It is hoped that the work of this 


Chesis Wil lwassist am promoring the practice Ger, this 





valuable management technique generally, but particularly 


Mobeni ume, Uo vo. Navy supp Ly sCOr ps. 
General 


It seems paradoxical that World War II which caused 
so much waste, both in human and economic resources, could 
also cause the development of new theories and techniques 
designed to indicate the most efficient and economic utili- 
Zation of resources. Historically, military necessity has, 
nowever, led not only to the development of new kinds of 
hardware, such as airplanes and digital computers, but also 
to the development of new scientific disciplines and con- 
ceptual understandings. 

One such outgrowth of World War II is the new discip- 
line known as "operations research" which has for its pur- 
pose optimizing the use of resources or, as stated by Jenny, 
"Operations research is concerned with optimizing the per- 
formance of a system. This requires the application of 
scientific methods, techniques, and tools." Within this 
concept, linear programming is a technique or tool of opera- 
Clone Trecscearcn. 

In the literature there are a great number of formal- 


tea. Gehinle1ons, and antormal descriptions), which attenper 


leoulding, RenneGiwine oc. len Io pve y «et cle 
Linear Programming and the Theory of the Firm. New York: 


Poe) Weacia il Lameco.g., USCC mip. oe. 
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© 
meme ossilty linear programing. —Simee simpli iacd definitions 
of a complex subject are seldom of much real assistance in 
the understanding of that subject, most available definitions 
fail to adequately describe the full significance of linear 
pueCtammng. “Among tne several classificarions to be Eound, 
linear programming has been variously described as ". . .one 
of the most important postwar developments in economic 


theory. . ae see es DramneGie Gl CONOMLC Tehneory. ae and 


(ewes SoiiaieneNacical Computing tecinicue. . poe Lie 
sense Gach description is correct but, at the same time, 
restricts the perspective in which linear programming should 
be viewed. For example, linear programming 1s used in econ— 
omic theory and received its first great impetus there, yet 
it is also used in the physical sciences as well; linear 
Programming iS not. just one mathematical computing technique 
but is, more properly, a mathematical approach (to problems) 
utilizing several computing techniques. Linear programming 
may be thought of as belonging to a larger class of problem 


solving approaches known as "mathematical programming." 





SSR aEHENO Robert, Paul A. Samuelson, and Robert M. 
Solow, Linear Programming and Economic Analysis. New York: 
MNeSteweridls Book -Co.7 Je. 1956, ps Vil 
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VoaINGen oe; 20s Ae Ooguceio cro ince Programming sand 
ee Tneorny. OF Games. London: Methuen & Com sles. 1 o6Gy 
oe ae 
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VeajCGej,, o, Readings in Linear Programming, News Yor<: 
John Wiley and Sons, Inc., 1956, De. V- 





Henderson and Schlaifer state that 
mathematical programming is not just an improved way 
of getting certain jobs done. It is in every sense 
anew way. It is new in the sense that double-entry 
bookkeeping was new in the Middle Ages, or that 
Meechanization in the Office was new earlier in the 
Century or Ehak alcomatiOnoins Ehe plant is new 
today. 

Programming problems, generally speaking, have to do 
with the optimum allocation of limited resources to sat-. 
isfy specific objectives. Specifically, programming prob- 
lems deal with situations where some resources, such as 
personnel, materials, equipment, and land are available, 
and are to be used in one or more ways to produce one or 
several results or products. There may be restrictions 
OW, Gena vallabr iri y rer lOnevOrvmore Of. Ehe Inpureresources. 
there may be restrictions on one or more of the expected 
output results or products; or there may be restrictions 
@n both inputs and outputs. Within the restrictions, how- 
ever, there are generally a vast number of possible or 
feasible solutions to the problem. From the great number 
of all possible solutions it is desired to find the one 
best or optimum solution, which will maximize or minimize 


(as required) some numerical quantity such as profit or 


cost. As Henderson and Schlaifer have said, 


Pender con Alexander and Robert Schlaifer, "“Mathe- 
matical Programming: Better Information for Better Decision 
Making", Harvard Business Review (May - June, 1954) p. 73. 


ee reuD Of .Limpeed resources Mist be "shared among a 
number of competing demands, and all decisions are 
‘interlocking’ because they all have to be made under 
that common set of fixed limits. 
It can be seen, then, that programming is intimately con- 
RPeGtooww? Giwoeel Silom making, rim Ene Managemen: Ol “resources. 

As pointed out earlier, linear programming belongs 
Memencn larger Olacci 1eCaclon co. Mathematical worogqrenming 
and deals with that class of programming problem for nes 
all relations among the variables are linear. These rela- 
Ehonsiips Must. be Llinecar Doth in the constraints and in the 
Objective function which is to be optimized. The linearity 
assumption is inherent in the discussion of linear program- 
mang enc Tes =sOlutcion techniques, “ihnis means, Simply thax 
as one variable changes, other variables have some pro- 
portional rates of chance. If the linearity assumption 
cannot realistically be made, or if the functions cannot be 
made linear by some suitable transformation of variables, 
then, linear programming cannot be applied. 

Until a more rigorous mathematical statement may be 
made, let the following description, in a more-or-less 
abstract context, be a sufficient definition for discussion 
purposes: Linear programming is the mathematical theory of 
the minimization or maximization of a linear function 


otpid., Dp. 744 


S) 
Subjectruo inear <eCnseralnes. ferceals with, tne 2nteraction 
of many non-negative variables subject to certain restrain- 


mg COnaLeELOnS . 


Historical Perspective 


Problems relating to maximization and minimization 
(extreme-value problems) have interested man ever since he 
first became aware of the concept. Euclid was concerned 
with finding the longest and the shortest straight lines 
@iat ecoould be drawn frond pointe .cO the circumference Of Va 
circle. Dido solved one of the first (although non-linear) 
Ineximizacion problems by Shreding “che: hide of an ox Eo 
bOum dene wperirery of Carthage. 

One of the most significant developments in forward- 
ing the understanding and solution of maximization-mini- 
mization problems was the establishment of the branch of 
mathematics Known as the "calculus of variations" which was 
first studied systematically by Euler and later by Lagrange. 

Linear programming problems are of a type Known as 
“constrained extreme-value problems" in which there is a 
Stipulation that the variables must be non-negative (a value 
of zero being allowed). Classical methods of the calculus 
of variations may be used to solve constrained extreme- 
value problems with no stipulations regarding non-negative 


variables. Classical methods, nowever, are not applicable 





when non-negative variables are required since they do not 
Gavewany assurance that the values ,olt those variables will 
moe EUrh OUt tO be negative. 

Very few, if any, disciplines have ever evolved 
independently, and linear programming is no exception. A 
considerable portion of the mathematical theory of linear 
programming is drawn from the theory of linear inequali-. 
ties and the theory of convex sets which were formulated 
during the past century. 

An immediate forerunner to the development of linear 
programming was the work done in 1928 by the late John 
von Neumann, formerly professor of mathematics, Princeton 
University, when he formulated the mathematical theory of 
games ("Zur Theorie der Gesellshaftspiele," Mathematische 
Biba len  §928 7 Vol. LOO) Dp. ASS EC(OOM IY Later, in 1944, this 
theory was expanded in collaboration with Oskar Morgenstern, 
Professor of Political Economy and Director, Economic 
Research Project, Princeton University, (Theory of Games and 
Economic Behavior, Princeton, New Jersey: Princeton Univer- 


sity Press, 1944.)° 


Riley, VeGasenC poduinwiie Gacs,, linea mb rode amint hog le 
Associated Technigues, Chevey Chase, Maryland: Operations 
Reccarch Office, The Johns HopkinssUniversity, 19585. .207. 


Dicer 2 OGE 





in 1941 Frank L. Hitehncoek, Professor Emeritus of 
Mathematics, Massachusetts Institute of Technology, pre- 


sented the original statement anda solution of the so-called 


"transportation problem", ("The Distribution of a Product 
maomtoeve talesOurces, TO Numerous Localities.” Journal oF 


Mathematics and Physics, Massachusetts Institute of 
Becnnology, August 1941, vols 20, no. 3, pp. PDA= 2 s0n 
The “transportation problem" (sometimes referred to as the 
"distribution problem") is one of several types of linear 
programming problems. : 

The second type of linear programming problem to be 
stated is what has been called the "diet problem" and was 
first presented in 1945 by George J. Stigler, Professor of 
Economics, Columbia University, ("The Cost of Subsistence, " 
JOurNelor Farm Economies, ‘May 1945, vol. 27, nes 2, 
pp. 303-314, )?° 

Each of these forerunners of the general linear pro- 
gramming model dealt with the optimization of a linear 
Runction subject to linear constraints but were not then 
specified to be linear programming problems. 

The formulation of the general linear programming 
model was, however, primarily the work of George B. Dantzig, 


Marshall K. Wood and their associates in 1947, while 





cee p. 284. 


1Orpid., p. 400. 
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emipleuwecd Oy 7ene Ue S. Aig Force. sec ehatetime tne Ul S.2aie 
Force was interested in developing a scientific basis for 
Mba lac eich BUCgGe tig ACC LSlONS 7 whe Tesule was the 
erganization of project SCOOP (Scientific Computation of 
Optimal Programs). Besides the immediate advantages accru- 
jig to the Air Force Erom the work of that research group, 
moro ject SCOOP contributed the “key” to linéar programming 
development and extension. The "key" was the development, 
by Dantzig, of a systematic procedure for solving the 
general linear programming problem; this procedure has been 
called the “simplex" method. 

The simplex method of solving linear programming 
problems gave the necessary impetus to research in this 
Field to allow linear programming to become a very impor- 
tant and valuable tool of modern theoretical and applied 
mathematics. Dantzig is quoted by Charnes and Cooper: 

The simplex procedure is a finite iterative method 
which deals with problems involving linear inqual- 
ities in a manner closely analogous to the solution 
of linear equations or matrix inversion by Gaussian 
Climinatiem.s «. « the term 7 simplex” evolved trom 
an early geometrical version in which (like in game _ 
theory) the variables were non-negative and summed 
eo unity.1 

Vajda, also quoted by Charnes and Cooper, states: 


tlonarnes, A. and W. W. Cooper, Management Models 


suc CD OuUStE ral Applications of Linear Programming, Volek 
New VYOnrr= * JON Wiley & Sons. Ine.) 196Gho) >. 242. 





de: 
The name (simplex method) derives from the accident 


that one of the first examples to be tackled by this 
method contained the constraint 


1 

which is the equation of a simplex in n-dimensional 

geometry. The name (simplex method) is now used for 

the procedure whatever the form of the (linear) 

constraints. 

The simplex method, as originally devised by Dantzig, 

has certain drawbacks, such as being very time consuming 
and having the tendency to accumulate round-off errors 


GB 
during the computations. What has come to be Known as 


the "revised simplex method", which is a streamlined ver- 
sion of the original procedure, was designed by DantZig, 
Orcharad—Hays and others at the RAND Corporation, in,order 
to overcome the disadvantages of the original simplex 
method, 

One other concept, that of duality, helped to pro- 
mote the growth of linear programming by indicating how cer- 
tain problems, otherwise intractable, could be solved using 
the theory of the Dual problem. It appears that Professor 
von Neumann first recognized that the dual problem existed, 
but it was Gale, Kuhn and Tucker who conducted the first 
investigations into duality and proved the first duality 


theorem. 





les nid. 
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Senremporary Perspect ve 


It became clear soon after the publication of the 
research report resulting from project SCOOP that the 
applications of linear programming far transcended only 
military application. Economists began, almost immedi- 
ately, to apply linear programming to Leontief's input- 
output models of economic systems; mathematicians began eS 
investigate the connections between the theory of games 
and linear programming; others have vigorously investigated 
wWarlous industrial applications of linear programming. 
opt readetOns Lange Erom tne-orlginal elassical applica— 
tions, to distribution and transportation systems, to the 
study of Kirchhoff laws of electrical-network cone 
and plastic-limit analysis in the design of Reese 

The rate of growth in the field of linear program- 
ming has been phenomenal and the theory is only in its 


infancy; today research is continuing at a tremendous rate, 


with no apparent end in sight. 
Reino) Ou mene Mucerna cure 
The very large quantity of published literature in 


the field of linear programming is almost overwhelming; 


oe DPe.-oce. 


Shp aay Dec O40. 
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for example, in 1957 Riley and Gass compiled a biblio- 
graphy*> SOntaimaing welsover 1one sehousania relLoreneces aid 
Since that time the available literature has increased con- 
siderably. Time has allowed for a review of only a small 
portion of the works available. 

it has been noted that, as in most fields of study, 
inconsistencies among the authorities appear in the liter- 
ature of linear programming. The most prevalent incon- 
Sistencies take the form of: (a) variations in the mathe- 
matical symbolism used, and (b) variations in the criteria 
used in solution methods. These inconsistencies are not 
serious and cause no real handicap; however, unless these 
differences are realized one who is unfamiliar with the 
literature may spend some uneasy moments in attempting ne 
correlate certain aspects of linear programming among the 
various authorities. 

In addition to the differences noted above, there is 
a wide range in the literature between the very theoretical 
approaches to linear programming and the “oversimplified" 
approaches; however, many authors attempt to arrive at some 
satisfactory degree of balance between the two extremes. 

PEEHOuUGM Luis dict roult: et Ccommene waniy One 
particular reference as better than any other, perhaps the 


a eens Semmens. 








is 
most concise and effective overall treatment, from both the 


theoretical and practical standpoints, are the books by 


Qe ee and Charnes and Cooper. / The book by Charnes 


and Cooper is more comprehensive than the one by Garvin and 
it contains some of the more esoteric approaches to, and 
moe catlOno OL, Linear programming. 

For anyone interested in a straight forward exposi- 
tion of linear programming and some of its industrial 
Boel iGaelOns ;Wwith a Manimum Of theoretical appendages the 
book by Ferguson and ear cene - is recommended, although it 
is advised that the lack of theoretical explanation may 
ultimately become something of a disadvantage to the 
meader . 

For those interested in a linear programming 


approach to economics the most celebrated and comprehen- 


SlVerWworm 2S that by Doriman, Samuelson and Solow.” In 


Sitcrsaie cacegory, DUC almGgst GServculy “entice tnecoreelcal 


level, is the book by eer ae 


Se evn Webter Veet orOjUeelon TO Linear yero- 
Seam a Neow Louk.” MeGray—ta li kood Co. incest o6U. 


1’ onarnes. Op. [Cut . 


Ser Gucony Robert O. and Lauren F. Sargent, Linear 


Peogmamtlng:, Pundamentals” and) Applications, New York: 
HeGreaw=nie ik BOO. COu., ne «956% 


Dorfman, Robert, Paul A. Samuelson and Robert M. 


Solow, Linear Programming and Economic Analysis, New York: 


MeGlaw Hii, BOok. Cow, nc. 1958. 


20 ae , ; 
Gale, David, The Theory of Linear Economic Models, 


mum eee ee Ee Se 


New wOrkKe. MeCrawenil) “Boek Co., anes, 960. 
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Several, rather short, treatments of linear pro- 


Gramming exist; probably the best are two books by 


Vajda.7*' ae Greenwald*> offers a short treatment of the 


Simplex algorithm, as does Ficken.** The presentation 


given by Greenwald is somewhat more readible and less 
theoretical than the work by Ficken which is almost 
exclusively a theoretical exposition of the simplex 


algorithm. 

The literature, represented in the periodicals, 
iS even more extensive than that represented in book 
form, and no attempt has been made to survey more than a 
very small fraction of articles available. 

Good, readable, general accounts of linear pro- 


Gr anming, including some Of tS sapolications may be found 


in the article by Henderson and Schlaifer;*> the article 


ence fu, Le hOCU Ge lO) eer Games. Opes Cites 


eva ida, REACRNGS wees (Oo Wohin OD wele. 


Sey eenwala, bakovea .Ulr ier, ioe ren Programming: 


Mee sD Lana ton Of the Simr lex Algor EieiaraG ad New York: The 
Ronald Press Gon OS 

Ceate on Eee ee SAD ex Metnod OF anedl rao] 
gramming, New York: Holt, Rinehart and Winston, Inc., 1961. 


*~ Henderson enceechl ate Oper ele. 
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by Cooper and Charnes; the article by Acrivos; and the 
article by Roconcers—- A good treatment of linear pro- 


gramming as it may be applied to long-range planning is 


Found in the article by Rapoport and Drews.°” 


eOceoner , Wiliam Wo. -ana-Abraneam= Charnes . "binear 


PaccueanwUiliag. Socleneciiic American, Auguste, [9547-vol. 19), 

Mon .2. 5 bo. 21-23. . 
ae Acrivos, “Linear Programming: How Does it 

Merce.) Chemical Engineering, ~6370) August, 1956), pp. 

pi o—Zlo. 

-eRosander, pee Greg ne. UseyOr iinears rogram inc 
HOoamprove the Quality: of Decisions,” Industrial. Ouality © 
Ser olw 2-9 etieren, W956. poe Ll=16- 

= RenopoLe, Leo A. and William P. Drews, "Mathe- 
matical Approach to Long-Range Planning," Harvard Business 
Review, (May - June, 1962), pp. 75-87. 
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THE LINEAR PROGRAMMING PROBLEM 


General 


An impressive range of management problems have been 
Bolved by standard linear programming techniques. Also, 

a great deal of research has been successfully FeconOLLsined 
to determine simplified methods of solution and reduction 
(by approximations), of some previously intractable pro- 
blems, to linear dimensions. 

The apparent basic simplicity of linear programming 
theory together with the wide diversity of successful 
linear programming applications may be deceiving to those 
not fully aware of the complexities involved. Recognizing 
linear programming problems to be existent in raw data, and 
translating actual physical situations into the language 
and form suitable for operations with linear programming 


techniques is not a simple matter. 


& 


Be COcna Zing Ulmecamarrogranning. Problems 


SLEUaLILONS Which Cone wWienin tne area Suitable for 
treatment by linear programming methods may be grouped into 
three broad categories, although there is only one abstract 
Statement of the Linear programming problem, as will be 


presented below. 





sey 

Wine SeatCOOrtes ©. Linear Ppreedramming (orOoblems are: 

Group tigi «lots “Gacegory 2nelvaes Situations, an) waicH 
the requirements for resources exceed the availability of 
the resources. Here, the problem is concerned with select- 
mag che particular requiremencts which can use the avail— 
able resources most effectively or economically. 

Cee UOw ie Mat Seeate deny weet sal Norco sl thae tors 
Weve Couci TeGuITements LOL ~ana avallability GL Tesourece sc. 
Problems of this type are solved if all demands are satis- 
fied and all resources are used; however, aalocaeion Gr 
the resources must be made in the most effective, or 
economical way when satisfying all the demands. 

Seely ie HOCH IGeOmlLi ee tle Category are wol celal ors 
where the availability of resources exceeds the demand for 
the resources. Problems of this nature require that the 
most economical assignment of demands, to certain available 
resources, be made. 

Re@arailess (Or any Of ene above three broad, general 
Categories into which a particular situation may be classi- 
Eleceseme ToOllowing aaqditional, factors Must also be tevidens 
in order to formulate a linear programming problem: 

(a) There must be some objective which can be clearly 
stated and may be optimized and which can be expressed, in 


mathematical terms, as some linear function. 


an Sane 
a. one beeen 
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(b) There must be restrictions on the extent of 
attainment of the objective and these restrictions must 
be expressible, in mathematical terms, as a system of. 
linear equations or linear inequalities. 

(c) There must be alternatives among which to 
choose, 

(d) The variables must be interrelated. 

Once a problem has been determined to be a linear 
mrocramming problem, in accordance with the above cri- 
teria, one may proceed to formulate the problem in mathe- 


Macical language - 


Problem 


It has been previously pointed out that linear pro- 
gramming deals with the interaction of many non-negative 
variables, which are subject to certain restraining con- 
ditions, for the purpose of determining maximum or minimum 
values. Some authorities state the standard linear pro- 
gramming problem as one of maximization (for example, see 
Saaty.°°)- Other authorities state the standard problem 
in terms of minimization (for example, see Churchman, 


30 : = 
seaaty, tomas ly) Matnenmacical Mecnods OF Opera— 


tions Research, New York: McGraw-Hill Book Co., iInc., OSS, 
Poe 26/—-166. 
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Ackoff, and Arnoff.~"). The latter generalization will be 


made here; then, the mathematical statement of the general 
linear programming problem may be stated as follows: . 


Paina Whe “values Of eee ee Wich 
1 Z 3 ia} 


minimize : 
n 
Pech Ge tee terete) ee ey = nee ee 
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subject to the conditions that: 


and 
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Churchman, C. West, Russell L. Ackoff, E. Leonard 
BeQOtG, 16 cla, Introduction tovOcerations Research, 


New TOsk: JVOnn Wiley and Sons, ince, 1957. op, 326-327; 





Zo 


OSs 


— 


given the column vec 


os, 





ye eacly 


weil). 


i = 2. -- 





Se A AR TR PES TT I ee a AE, ON RE EA A OS Ne a ay ey iv 


a a ® ® ) ea | 


el 








Zi 


Getermine values of xy. Koy Kay. eee x which minimize the 


iinear euneccione is 

n 

3 CX, = C1 X)+C5X5tC3% 3+ ee ce ae — A OG oy ue een On 
j=l 


sub ject to the conditions that: 


j 
n | 
oF a ae = PX) TPHaX5tP AMT ~ - + P tn = ie : Sy wit pm) 
teat 
Where : 
X, are unknown variables, and (4 = 1,2,3,22-,M) 
a.., b., and c, are given constants (5 = ie ee 
1 J 1 4 


Mhe proplem 1S to select, out Of the inkinite number 
Sa rcOlueilons Of Gquaticons (2.lce) or (2. le). the set or solu 
tions which contains only non-negative variables and for 
Wihemwene “linear form of equation (2.,la) is a ‘minimum, fn 
other words the central problem is to minimize equation 
(2,.la) subject to equation (2.1b) and to the indeterminate 
System of equations (Z2.1c) or (2.ie). Equations (2.lc) 
are the constraints of che system, equation (2.1b) represents 
the non-negative conditions and equation (2.la) is the 
Objective function. 
€ 

Thus the general linear programming problem has been 


Stated in its standard (or canonical) form, where all the 





Ze 
constraints are equations, all the variables are required to 
be non-negative and the objective function is to be mini- 
mized. 

Linear programming may also apply to situations 
where the constraints, instead of all being equations, are 
all inequalities or a mixture of equations and inequalities, 
or some variables may be negative, or the linear functional 
is to be maximized instead of minimized. Any of these 
alternate situations may be reduced to the standard form of 
linear programming problem as indicated above. 

When a constraint is stated in the form of an inequal- 
ity it may be changed into the form of an equation by adding 
a non-negative slack variable to the left hand side of the 
equation. For example, if we are given the less-—-than 


inequality: 


by adding a non-negative slack variable (S.) toe che nlvert 
hand side of the inequation it is transformed into the 


following equation: 


Similarily, if we are given a greater-than inequality: 


> 
Sool ee ee i te 
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the inequation may be transformed into an equation by adding 
a non-negative slack vector (S, ) to the right hand side of 
the inequality (or, transposed to the left hand side it 


becomes [-S,] ): 


If some variable is not restricted to non-nega- 
tivity (i.e., an unrestricted variable) it may be trans- 
formed into a non-negative constraint to suit the standard 
linear programming form. Since any number may be expressed 
as the difference of two non-negative numbers the variable 


can be defined in terms of two new variables: 
> ae tee? ae — ae Gar ao, See ee C2 2. 


i ii 
wheres: Ba Oana _ 2 O 


Substituting the above relationship (2.2) into equation 
(2.1c) or (2.le) the system again conforms to the standard 
cOrin. Hadley~* offers an interesting demonstration that the 
value of Xx, uniquely determines the values of X and 4 


lf, instead of wishing to minimize, it is desired to 


maximize an objective function so that: 





eereaiiey, Go, uinear Programming, Reading, Mass: 


Addison-Wesley Publishing Co., Inc., 1962. p. 169. 
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Men, tniis relationship may be rewritten: 


= |p} c 4X, = -c)X, ~ C5X, - CX. - 2... cox, = 54 via) 


wa jda states: 
(Wie aEOoLeMmeon MAXIMNLZIing an Objective. TunctLon is] noe 
fundamentally different from that of minimizing one, 
because maximizing is equivalent to minimizing the 
hegazive, 

Thus, the maximization problem is transformed into a 

minimization problem of the standard form. This procedure, 


however, is not identical to the duality concept discussed 


below. 
The Mathematical Model 


The concept of a mathematical model is intimately 
associated with the subject of problem formulation. In the 
preceding section we have described a mathematical model 
of the general linear programming problem in its most 
abstract sense, 

A mathematical model is a device which aspires to 
be a replica of a system which is to be studied in some 


detail; without constructing a mathematical model one cannot 
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formulate the linear programming problem. Mathematical 
models may vary in the degree of faithfulness with which 
they adhere to an actual physical situation, depending 
upon the desired accuracy which is expected or the pur- 
pose for which the model is intended. 

Once a mathematical model is established it may be 
studied from various perspectives, utilizing various 
assumptions in order to gain insight into the “real life" 
circumstances of the problem. The versatility of a mathe- 
matical model will be demonstrated in Chapter V and Chap- 


ter VI. 


Mesa imal Proplem 


It is convenient, when discussing linear programming 
problems, to speak of the "primal" problem and its “dual" 
problem. 

The standard problem statement as indicated above 
1s generally considered the statement for the primal prob- 
lem, although it is really immaterial which problem state- 
ment is called the primal and which is called the qual. 

To facilitate discussion of the primal problem 
and the concepts presented thus far, some simple (bi- 
variable) examples will be used for illustration. Examples 


with only two variables will be used, obviously, in order 





to allow for graphical representation, since bi-variable 


Petacionships may be presented on cartegian coordinates. 
(Example Problem No. 1) 


Peacsiiomeicde a WU. oo Navy Suoply Ofricer* inds, thac 
he must order a quantity of gasoline for stock and for 
immediate use. He has available storage capacity for 
?,600 gallons. He normally issues two grades of gasoline 
in about equal quantities (each grade of gasoline has 
approximately equivalent demand rates). The supply officer 
Knows that he has an immediate requirement for at least 
in coo Gallons Of either kind, or a combination of both 
Panosmor CasOline,. but due to special operational require= 
ments all of one type of gasoline (x, ) plus one-half of 
the other type (X,) must equal or exceed 800 gallons. The 
officer has been informed that there are only 700 gallons 


of type X5 enc 300 gallons OL tyoe xX, available co fill cme 


th 
Order at the present time. The gasoline costs 15¢ per 
gallon for type x, and 13¢ per gallon for type Xo - The 


supply officer must determine what quantities of each type 
of fuel to obtain in order to minimize the cost to the 
Sever ument . 

In order to formulate the problem in terms suitable 
for a linear programming approach, a mathematical model has 


to be constructed, as follows: 
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Me "Objective woicn ls Gdesired is Ene Minimization 


ee COST, Or: 


Oe aXK 
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moitecn conforms to equation (2ila) of the standard linear 


programming problem. 


Since negative quantities of gasoline would not be 


a reasonable expectation, 


wer may State: 


IV 


(Ze) 


NV 


These relationships (2.3b) correspond to equation (2.1b) of 


the general linear programming problem and also indicate 


EAat a Graphical representation of the problem in cartesian 


coordinates will be applicable in only the first quadrant. 


fier condi ELGns Ob -COnStraints of Jtne problem, 


Sorrescponaing to the system of equations represented by 


relationship (2.lc) of the standard linear programming 


problem, may be 


1 2 
Ay ar Ao 
xy - 3K 5 
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stated mathematically as follows: 


IV 


IV 


tA 


A 


1,600 


1,000 


800 


S00 


700 


CSGoreage ICa ogee 7 ) nee vluves ce aa 
POrder Sousa it 7). Mo.face Zee Sie 
(Special Usage Requirement). . 
Pipaim CaciOn Of % Availability) 


[ Limptat)om OF xX, Availability) 


(2356) 


2 Sco) 


(2.3e) 


Score) 


(220g) 





28 

Utilizing figures (2.la) through (2.le) we may 
demonstrate the meaning of the constraining inequations; 
in figure (2.1f) the meaning of the objective function will 
be shown. 

Equation (2.3c) describes the storage limitations 
in terms of a less-than inequality. This condition is 
indicated by the line in figure (2.la) at the maximum 
limit (extreme); the arrow shows that this constraint may 
take on an infinite number of other positions toward the 
Srigin of the graph. 

Equation (2.3d) describes the minimum order quantity 
and is expressed as a greater-than inequality. This 
relationship is shown by the line in figure (2.1b) at its 
minimum limit (extreme), with the arrow indicating that 
this constraint may take on an infinite number of positions 
away from the origin of the graph. 

Scimplaiiy eGuations (2.36), (2.36) and: (2.30) ane 
shown by lines at the extreme limits of the constraints in 
figures (2.1lc), (2.1d) and (2.le) respectively. Again, 
arrows indicate the direction in which the constraint may 
move to assume an infinite number of other positions. 

Figure (2.1f) shows the objective function at several 
of the possible positions out of the infinite number of 


Other positions which it could assume. Since the objective 
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So 
function is not constrained (except by the non-negativity 
of the variables) it may take on any position from zero, 
@eche Origin, out to infinity. 

Superimposing all the individual graphs one above 
the other we can see a consolidated representation of the 
problem, as a whole, in figures (2.2a) and 2.2b). The 
enclosed polygon in figure (2.2a) indicates the area of 
feasible solutions (shaded area); any combination of 
X) and X54 within the shaded area, would satisfy the con- 
straining conditions defined in equations (2.3c) through 
(2.3g). However, since we are seeking an optimum solution, 
the solution will occur at one of the extreme points of 
the convex polygon formed by the constraining boundaries. 

Since the problem is to minimize the objective 
function, imagine the lines (with a slope of -15/13) as 
shown in figure (2.1f), progressing from the origin out- 
MoiGdeaimwrvguire (2.20). The first point of tangency of 
this line with the polygon, shown in figure (2.2b) by the 
shaded area, is at the point X 


= 600, X, = 400. This 


ll 


point, then, is the maximum order quantity for minimum 


2 


expenditure. The cost of this order may easily be computed: 
(20 25 (G00) (GO 213) 400) 9 $1.42.00 
To make certain that this actually is the minimum cost we 


May examine the cost of the other extreme points. If the 
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c= 450, Xo =/ OU. ene Coste would be $159550. (1s 


the order were x4 = )510).0)" Xo = 200 the cost would be $146.00. 


order were X 


If it was decided to order the maximum of available 
quantities xy = 800, x5 = 700, the cost would be $211.00. 
Therefore, by means of the method first used it is seen 


that we obtained, directly, the minimum cost of the order 


which was the objective originally sought. 
(Example Problem No. 1 [Modification No. 1]) 


Continuing with the example used above, assume that 
the supply officer had instituted a policy that there may 
be an imbalance of stock in favor orf larger quantities of 


Xo, Since it is somewhat more difficult to obtain and it 


is less expensive than X The officer, however, has 


1° 


Seeciried that x5 should not exceed x) by more than 200 


gallons. The new conditions, required by the policy, may 


be shown as follows: 


eG Ae (Minimum imbalance of stock). .(2.3h) 
ee X5 Z10) 
X5 - xX = 5760 (Maximum imbalance of stock). .(2.3i) 


= tx = 00 





S. 
These additional constraints, along with the original con- 
straining conditions are shown in figure (2.3). It can be 
seen that a new polygon of feasible solutions is formed, 
Proceeding as before, imagine the same line, representing 
1+ 0.13X, = Z) to be pro- 


meessing trom che Origin Outward. The first point of 


the objective function (0.15x 


tangency is now: x1 = 2ooeoo x. = 533.33. The new 
minimum cost of the order, under the constraints, is 
p49 .33. 

b One would imagine that the supply officer, by 
instituting such a policy, as described above, was attempt- 
mig CO Maintain erlticient economical operations. it is 
readily apparent, though, that he has been successful only 
in equalizing his order for the two types of gasoline 

while increasing the cost of the total order, 

Maietlartcer example points @ut one Of the many, bus 
probably least obvious, advantages of linear programming: 
the evaluation of policy. When certain criteria are con- 
Sidered essential to operations these criteria can be 
evaluated in terms of their additional cost, as shown above. 
Perhaps, when it is shown exactly wnat certain policies 
cost, 1t may be determined that their "essentiality" is 
only imaginary. 

Linear programming has found many useful applications 


in private enterprise. Here, again, as in the example given 
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pmove, One May be seeking tO minimize costs; however, unlike 
Navy or governmental operations generally, which are uncon- 
cerned with profit “per se", one most probably would be 
more concerned in private industry with maximizing profits 


of the business, 


(Example Problem No. 2) 


Consider the following very simple example: an oil 
company manufactures two types of lubricating oil from one 
fomna OF crude Oil stock. Each gallon of product requires 
equal quantities of one gallon of the crude oil, but the 
crude oil is available only to the extent of 12,000 gal- 
fons per month and it costs 5¢ per gallon. In order to 
make the two products a certain machine processing opera- 
tion is required where the cost of operation is $2.00 per 
four. Product X, requires one-half of an hour of machine 


a 


time operation per each gallon produced; product Xo 


requires one hour of machine time to produce one gallon of 
product. The machines are available for production of these 
ewO products for only 500 hours of each month. If the 


company can make a profit of 54¢ per gallon on product xX, 


and 75¢ per gallon on product X., what amounts of these two 


Z 


products should be produced in order for the company to 


receive the maximum profit? 





3c 
First, the problem may be stated mathematically in 
consistent units (i.e., monetary units in this case): 
crude oil available: (12,000)(5¢) = $600.00 
machine time available: (500)($2.00) = $1,000.00 


X, machine requirement coefficient: (4)($2.00)=$1.00 


X., machine requirement coefficient: (1)($2.00)=$2.00 


The objective function for maximum profit is: 


iu 0.075X. ah aA re edo i tae Mae. eee ees) 


Negative quantities of the products cannot be made; 


O2055xK 


therefore: 


rs 
NV 
oO 


ee ee ae ee eg 8!) 


IV 
Oo 


mane constraining conditions are: 
hoe = 600 (crude oil limitation) . . . .(2.4c) 


x1 + ZX =A, 000 (machine processing limitation(2.4d) 


The convex polygon of all feasible solutions is shown 
by the shaded area in figure (2.4) which graphically 
represents the above relationships. Since the objective 
in this case is the maximization of profits we must imagine 
that the line representing the objective function is pro- 
gressing from infinity toward the origin (at a slope of 
~0.055/0.075). The first point of tangency to the convex 
polygon, which is the maximum DOLnt, 2S ae = ZOO; 


1 
Ko = 400. The maximum profit is, therefore, $41.00. 





ew 
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It is interesting to note row the various profit 
figures affect the Ob jective Baik On wa hOenGw ot Le eee 
determines the optimum quantities o*7 production. 

Figures (2.5a) and (2.5b) reproduce the convex polygon 
Trormed by the constraining relationships of the previous 
example and with various new objective functions having 
Pea Gea serOri.: Margins. 2 facie (2.54) the solid line 
EP PrescomecrvanvOpjeceive fUunck2O94O. equal protiecabi lity 
between the products X, and X. (:or example, assume: 
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represents a profitability of xX) slightly exceeding xX, 


(for example, assume: 0.0605X, + 0.060X., = Ze Note thas 


2 
in the case of equal profitability the line representing 

the objective function would come to rest directly upon 

the line representing the constraining relationship describ- 


ing the limit of availability of the crude oil (X) + Ke SoOC je 


Z 
This situation indicates that any feasible solution on that 
line (including the maximum solution of the previous 


| * 
example) would be equally profitable. Possible solutions 


in this instance are: 


X, = 200; X, = 400; Z = $36.00 
X, = 300; X, = 300; Z ='$36.00 
X, = 400; X, = 200; 2 = $36.00 
X, = 500; X, = 100; Z = $36.00 
X, = 600; X, = 0; Z = $36.00 


Determining such a fact as this might have important bearing 
on management decisions since all elements necessary fOr 
final effective decisions may not enter into any basis of 
the problem solution because of certain non-quantifiable 
factors concerning the product data. For example, it might 


be more desirable, for good labor relations, to produce one 





* 
Note that intermediate solutions such as: 


ee Oe lve 


values are acceptable. 


aos 349.9, are equally valid if non-integral 
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agreeable product rather than some other product which may 
be offensive to personnel for some reason; or, one product 
might be less hazardous to make than another product, 
Knowing that equivalent profitability could be gotten from 
various product mixes, management would be able to make 
more valuable decisions. Under certain circumstances, 
however, aS in this example, when a slight fluctuation in 
the profit margin occurs (like the incremental change shown 
by the dashed line of figure [2.5a]) there is only one 
optimum product mix. 

Figure (2.5b) represents substantially the same idea 
ee equals prOritability but with regard to a profkit margin 
differential between the products. The solid line repre- 


sents the relationship of product X. being twice as profit- 


Z 


able as product X, (for example, assume an objective 


a 


PmiMCeron~=— OLSO5X -+ 0, 10K. = 2). Here the objective func= 


Ht 2 


tion would come to rest upon the line representing the 
constraining relationship resulting from limited machine 
capacity. Again, we have a situation where we may have an 
infinite number of solutions (assuming that the solution 


may contain non-integer results) anywhere from eae 


xX -OUUPUP te and including X-~=.200<- xX. *= 400. in any. 


al 2 


case the profit return would be a constant $50.00. 
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Next, if we allow a small increment in the profit 
mergin Of product X5 (assume, now, an objective function: 


emgox. + O.,101X, = Z). In this case the dashed line, 


Li 2 


Representing the objective Function, first touches the 


polygon of feasible solutions at Xy =r X5 = 500 indicat= 


mie that maximum profitability can be made Erom production 


mimmonly oroduct X. (with a profit return of $50.10). 


Z 


iveene same manner, had product. X- “increased slightly 


ie 


in profitability instead of product X, (assume an objective 


Z 


1+ 0.10x., =~2)20 The maximum) prot ie 


emen would be $50.02 and occur only with the production of 


mamction Of: O.O501X 


xy =1Z00l- X5 = 400° and this was the original solution 
obtained. 

In addition to showing how various product combina- 
tions may, under certain circumstances, return identical 
profits the foregoing illustrations indicate how some slight 
Merlo @n i) ProlLit Margin May become very critical-to the 
decision making process. This points out that seldom will 
a solution to a given problem become static. Continual 
vigilence and review of pertinent conditions is obviously 
necessary to insure that optimum solutions, once gotten, 


Seman Optimal). 
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The Dual Problem 


To every linear programming problem involving maxi- 
mization there exists a specific corresponding minimiza- 
tion problem involving the same data and having the same 
Optimal solution. The first problem explained above, is 
called the primal; its counterpart is Known as the dual. 

As mentioned previously, it is really immaterial 
which problem is considered the primal and which problem 
is considered the dual. It is equally valid that for every 
linear programming problem requiring minimization there 
also exists a specific maximization problem involving the 
same data and having the same optimal solution. 

Utilizing the data of the previous example of profit 
maximization, considered in the last section, the pertinent 
data is displayed in table (2.1). It can be seen from 
table (2.1) that the primal problem is formulated by taking 
the data row by row whereas the dual problem is formulated 


by taking the data column by column. 
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THE PRIMAL/DUAL PROBLEM 
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SOLUTION 
$41.00 





PRIMAL PROBLEM DUAL PROBLEM 
Cy X, + CX, = Z (maximum) b,W, + boW, = 2% (Manion) 
a X,. + aX = a..W, + a, W = 
ih aes. eee = bes ee eal 
a » an or. eee 4 ie a..W. + awAW Z 
a 2D 9 ein DO ee 


Specifically, the primal and qual problems, for the given 


example, are formulated as follows: 


PRIMAL PROBLEM DUAL PROBLEM 
0.055x, + 0.075X., fe uice es) SOOW, + 1, 000W., = 2% Mal) 
wees Ce Solon se W. + W. = 0.055 
1 Dia ; 1 7 ea 
< > 
X, + 2x. = 1,000.00 Ww, + 2w, = 0.075 
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Graphically, the primal problem was presented in 
figure (2.4). The graphical representation of the dual 
problem is shown in figure (2.6). It can be seen in fig- 
ure (2.6) that the polycon of feasible solutions extends 
from the lines representing the constraining relation- 
ships out toward infinity. 

In the dual problem the objective being sought is. 
to find the minimum value and, as before, imagine the line 
mepresenting the objective Eunction to proceed from the 
origin (with a slope of -600/1,000 = -3/5) toward infinity. 
In figure (2.6) it can be seen that the first point of 
tangency of the line with the polygon of feasible solu- 
Selous 1S at Ws = tO 2055 Ww. = U5 O22r * tne bya tue Om cic 
PojeGtive function at this point is: (600)(0.035) + 
(1,000)(0.02) = $41.00 which is the same solution reached 
ki the maximization of the primal problem. 

The proof of the duality theorem, which is somewhat 


involved, is given by Ces 


Conuentson Frrooliem Formulation Context 


ELL SS Ay ES SSS EE 


Pee snoulra be “pelited out “ab this poinG®, 1Li-it 2s 
not already obvious, that the context of a particular prob- 


lem is really irrelevent (except in the matter of final 
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GRAPHICAL SOLUTION OF THE DUAL OF EXAMPLE PROBLEM NO, 2 


interpretation of the results of the solution). What is 
important and indispensable is the proper construction of 
the mathematical model in the beginning of problem formu- 
lation. When a problem is stripped of its verbiage anda 
mathematical model is properly constructed in terms of 
interrelated variables (say: Kir Xone 2 2, x) it does not 
matter what the variables really stand for (assuming they 
have been properly defined at che outset). For example, 

in the foregoing problem, where profit from the manufac- 


ture of two products was maximized, the problem could have 





4 
been stated, just as well, as follows: AU. S. Navy Supply 
Officer has a disbursing clerk who can either compute pay 
or audit vouchers, or do both jobs sequentially. The clerk 
can process one pay record (on the average) in 3.30 min- 
Mees- tO audit one voucher it takes the clerk 4% minutes 
(on the average). Both operations require the use of a 
Mmepe-woiter- at the rate Of One minute per document fEor both 
pay records and vouchers; there are 10 hours of typewriter 
time available to the clerk during one week. Each opera- 
tion also requires the use of a desk calculator at the 
respective rates of one minute per pay record and 2 minutes 
per voucher; a total of 16 2/3 hours of time are available 
to the clerk for using the calculator during one week. The 
supply officer desires that the disbursing clerk produce as 
much work as possible during the week, therefore, the 
officer must determine how to best assign the work to the 
clerk. 

In essence the problem just described is identical 
to the preceding problem although the context is entirely 
different; by constructing the mathematical model this 


becomes evident. Let: 


x4 pay records processed 


Xx 


5 vouchers processed 
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Since it is desired to maximize the clerk's weekly output, 
then the objective function is: 


Seo Cn tau OO KX. = FauiMax) 


dk Z 


Seoject to the restrictions: 


1/60 X, + 1/60 x, = 10 
1/60 X, + 2/60 x, 216 272 
or, 
0.055X, + 0.075X, = Z (Max) 
X, + X,, Seco 
X, + 2X, = 1,000 
ana 
Xx. = 0 
eS 2a 


Thus, we have the same problem as before but in an entirely 
different setting. In this latter oroblem, of course, the 
solution is exactly the same as before (X, oO. x5 = 400; 
Z= 41). The interpretation, from the model to the physical 
Situation is the only difference; the numerical data and the 


mathematical model are identical. 
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BASIC METHODS OF PROBLEM SOLUTION 


EMPLOYED IN LINEAR PROGRAMMING 


General 

Until the development of the simplex method, there 
were no general solution techniques available with which. 
to solve linear programming problems. Ferguson and 
Sargent report that “the simplex method is the basic linear 
programming method from which the other known methods have 
been serves Gale states that *he simplex method is 
" . .one of the most celebrated computational procedures 
in recent mathematical history.", and “. . «so successful 
has the simplex method been that from the literature one 
fgnie her lea to believe that linear -~programming 1s che 
simplex method, ""° 

Many modifications and extensions have been developed 
since the original simpiex method was first reported. A 


“revised simplex" procedure has been established and is the 


algorithm generally used in computer codes, but Hadley 





3° Perguson Glee ei Milnes Oe oO cle 


38cale, Ope Citi, De. 205: 








3) Ak 
Pepocca weiace che revisec simples method Nas not mec with 
wide acceptance for hand Soetoro | 

Another outgrowth of the simplex method has been the 
development of the "transportation method" (or, sometimes 
called the "distribution method"), end the "modified dis- 
Breibution method" (or, “MODI method"). 

The SuUserOr any Of Lhe sineaer programming, sOLUtiCn - 
methods requires the manipulation of a great deal of 
arithmetic, and the real power and usefulness of linear 
programming lies in the solution of very large mathematical 
models where the employment of hand computations would be 
marealistic, and only digital computation machinery can 
Seecucro!y perform the necessary arithmetic. There is a 
place, of course, for hand computations, when relatively 
small problems are being considered, and electronic com- 
putation equipment is not available. 

In the following sections of this chapter a dis- 
cussion of the simplex method and the modified distribution 


(MODI) method will be presented. 


Simolex Method (Minimizaticn) 


The literature abounds in examples of using the 
Simplex method to solve maximization linear programming 


problems. There are not so many examples of the use of 





ee alleys Oe, ee, ©, 2illae 








a2 
this method to solve minimization problems, however. The 
method of solution is essentially the same whether one is 
Bee, oesing Maximization Or Minimization problems; but there 
Mee Slight dirterences which ere not altogether obvious, 

In this discussion a minimization problem will be 
treated by converting it to a maximization problem and pro- 
ceeding with the simplex computations. Where there are 
differences in proceedure between minimization and maximi- 
gation, the differences will be pointed out in the ensuing 
discussion, 

The discussion of solution methods may be facilitated 
by the use of an example; in this discussion example problem 
No. 1 (modification No. 1), presented and solved graphically 
in chapter II, will be used. The mathematical model (equa- 
miONSsii2.sa |) through | 2.31]), introduced in Chapter 11, 


Ls reproduced below: 


Objective Function: 0.15x, - 0.13x, =o. ee a) 
Sceiseraints: xy + x, = imo O8 
x, + x = 1,000 
x, + 5X, = 800 
xX, =  g00 
X, =P 2700 
- X, + X, = 0 
- X, + X, = 200 
X. = O 





DS 

Prior tO attempcing a Solution of a system of linear 
inequalities (as shown above) the inequations must be con- 
mercea to equations as indicated in chapter II. This con— 
version may be accomplished by the introduction of slack 
variables. These variables accomplish just what their name 
implies: slack variables take up (or, are assigned the 
quantities of) any difference, or slack between the total 
Biantity, called for by the equation, and the total quan- 
elty assigned to the primary variables. | 

An inequality of the "less-than" ( =) form may be 
eonverted to an equation by the addition of a slack vari- 
able to the left hand side of the inequation. An inequality 
of the “greater-than" ( = ) form may be converted to an 
Povalieon py the eddition of a slack variable to the right 
hand side of the inequation; then, this slack variable is 
subtracted when transposed to the left hand side of the 
inequation. 

Since the slack variables merely absorb unused 
fe -OURGCS Ss ehey add Nothing to.the Brorle-and incur nothing 
in the way of costs. Therefore, slack variables must carry 
Beno © )MproOLl: and zero (0) Cost GCoctEicients> however, 
these variables may be assigned penalty coefficients under 


Cebteifecirecumseance s. 
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The mathematical model, shown above, consists of the 


following system, after the inclusion of slack variables: 


Sajective Functions 


O.15X, + 0.13X, + OX, + OX, + OX, + OX, + OX, + OX, + OX, = 2, {Mateo 
xX, + X, + X, = le OOd 
X, + ee - X, = 1,000 
X, + 2X5 = = 800 
x) + X- = 800 
X5 “- xX == 1OC 
~ Xx, + Xo - Ko = 0 
- X,+ xe + Xg = 200 


In all solution methods employed in linear program- 
ming one always proceeds toward a final optimal solution by 
starting with a feasible solution (which is not necessarily 
optimal but is a possible solution to the system of equa- 
tions). Generally, in the simplex method, a first feasible 
solution is gotten by allowing the structural variables 
(X) and x, in the present example) to equal zero. This is 
equivalent to saying that nothing will be produced and no 
resources will be utilized. 


Now, notice in this example that if we take the 


approach of allowing xX) and xX, to equal zero we cannot 





52 
obtain a feasible solution. This misfortune comes about 
because we would be saying, among other things, that: 

x 


4 


20 


=a OOG 


800 


and we have, therefore, violated the non-negative variable 
requirements (X = Ol: 

To overcome this difficuity we may adopt a device 
solely for the purpose of obtaining a first feasible 
solution. The device which will be used is the introduc- 
Lon OL “artificial variables” which have no physical rela- 
tion to our problem and, therefore, must not enter (at a 
positive level) into the basis of the final optimal 
BeluttOmen, TO -prevent the artificial variables Erom becom— 
ing a part of the final optimal solution we will enter them 
at an arbitrarily high cost (M) since we are attempting 
minimization of the objective function. Since the arti- 
ficial variables have no bearing on the problem and must 
not enter the final optimal solution the high cost (M) 
may be carried throughout the computations without being 
actually specified numerically, if it is always remembered 
that it (M) is the largest number in the system. 

The mathematical model, after introducing the arti- 


ficial variables, consists of the following system: 
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Ob jective Function: 


15x, + ;. 13X,+OX,+0X,+OX_+OX-_ +OX+0X, +OXg+MX, (+MX, 1 +MX, 5 = (Min) 
X,+ X,+ X3 =] , 600 
Xj + Xo - Xy + X19 i 91319, 
X,+ 2X5 - X- + Xy = 2cC 
X) + X, = ‘800 
X, + Xo =e Oe 
- X,+ xX. - Xo + X15= 0 
a X,+ Xo + Xo = . 200 


At this point we encounter the first difference 
between minimization and maximization problems. Above we 
have entered the cost of the artificial variables as a 
positive high cost (M). In a maximization problem we would 
have to enter the artificial variables at an arbitrary high 
negative cost (-M) in order to prevent their entry into the 
minal optimal solution. 

: In chapter II it was stated that “maximizing is 
equivalent to minimizing the negative," Hadley offers a proof 


GE thiss relationship. *° Numerically this concept may be 


demonstrated as follows: 





Sete sO Lome 





Si 


Biven the following series: 


then five (5) is the maximum. Negating the series we have: 
Pee =D 1, 0, HL, #2 
then, minus five (-5) is the minimum. 
Using this fact, as demonstrated above, we restate 
the objective function as follows: 


-.15X, -.13X.,-OX.,-OX , -OX, -OX, -OX..-OX,, -OX, -MX, ,-MX, , -MX 9=2 (Max) 


1 2 3 4 > 6 a 8 2 10 ileal 1 


The tableau format is the universal device used for 
linear programming computations and the present system of 
equations is shown in tableau (3.la); it is a representa- 


Meier neowarrste reasible solution (with xX .=xX2=X -=X_-=%2=Oue 


lia ae ie os es 3: 
Column (c,) contains the cost coefficients (c 5) of the 
respective variables in the current basis (these variables 
are found in the column headed: "vectors in basis"). The 
Main body of the tableau contains the detached coefficients 
ey” of the variables (where no entry is made, a zero is 
to be understood). As a convenience all the cost coeffic- 
ients of the objective function have been entered above the 
tableau; and the columns headed "Key Row" and "Ratio", as 
well as the row headed "Key Column" have been added for 
discussion purposes. 


Column (b), Row (2) represents 2 c_b: in the first 


b 


tableau this is computed as follows: 
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(0)(1,600) = 0 
(-M)(1,000) = -1,000M 


(-M)( 800) = 800M 


(O)}{ 800) = 0 
ej 700) =.0 
Ca) 0) = 0 
(Ope (200) 4=,6 


RO Geark -1,800M 


Sener values in row 2a) tor “tney veri ous -columc (X,) are 

computed in the same manner as above (i.e.,: oz Cay 3). 
ROw (2. - C5) is computed by taking the (z 5) as 

computed above and subtracting the respective column vector 


= i=—,1 > eenenwese = c 


cost Ser 5 for example, 2, = -M, c, l 1 
—-M - (-.15) = -M + .15. Normally row (Z 4) is not recorded 
since our main concern will be with row (2, _ c,). 


After an initial tableau has been constructed one.is 
ready to begin iterations, in a step-by-step procedure, 
Eewarad a stinal optimal solution. 

Step 1. To begin an iteration find the most nega- 
tive (2 _ Cie LdeNnti Ey Les: column (X 5) as the “key column" 
(note the double asterisk in the "key column" row of the 
tableaux). The x. will be the variable entering the basis 
in the new tableau, 


SEeP ico Divade Cach mumber in the "key column, <inro 


its corresponding value in column (bd). (Note “ratio” column. ) 





ae, 
® 
Find the quotient having the smallest positive value; 
identify that row (X.), associated with the smallest quo- 
tient, as the "key row". (Note the double asterisk in the 
"key row" column of the tableaux.) The variable (X,) is 
the variable which will be eliminated from the basis. 

In attempting to determine the smallest positive 
valued quotient it may happen that there are "ties" where 
two quotients have the same smallest positive value. In 
the case of ties the following procedure will allow one 
to continue the general procedure: 

a. If the "ties" are not caused by zeros in the 
numerators of the ratios then the ratio which has the 
largest denominator is chosen and that row is identified 
as the "Key row". One would then go on to step 3. 

b. If the "ties" are caused by zeros in the num- 
erator of the ratio, then the denominators of the ratios 
are each divided into each value in their respective rows 
Scateing ac the constant column and moving to the right 
comparing at each column. When one of the new ratios 
becomes smaller than the other, the tie is considered 
broken and the row containing the smaller value (i.e., 
smaller ratio) is designated the "key row" and one then 


would proceed to step 3. 
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Steps Conley ene nuUMber at tie Inter seqklonuor 
the “key row" and the “key column" as the "key number". 
(Note the circled numbers in the tableaux.) 

Such wae pivice aii numbers 1n- che = Key row! by che 
"key number" starting with column (b) and going through 
column (X_). These quotients are entered into a new 
tableau in the corresponding row, and that row, under the 
column heading "vectors in basis", is designated by the “ 
of the "key column" of the previous tableau; this is the 
entering variable. The cost of the entering variable is 
placed under column (c)) in its respective row. The 


memainaer Of rows under columns (c and “veetors in basis" 


»? 
are reproduced in the new tableau from the previous tableau. 
Step 5. The values for the remaining boxes (or 


cells) of the new tableau (including rows Lz.) and 


Tae - oe) as well as column [b!) are computed as follows: 


= corresponding corresponding 
value in Va lie: in “key row" x “key column" 
new =} old to $0 UNS Ey NuMper 








tableau tableau 
("key number" ) 


As an example, column X, of the second tableau (i.e., 


8 
tableau 3,2b) is computed below: 
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CORRESPONDING CORRESPONDING 


aa IN “KEY ROW" x "KEY COLUMN" _ ecet 
ROW NUMBER NUMBER = zd 
TABLEAU oS TABLEAU 
"KEY NUMBER" 
ee ae Z 
XK O - 7 as i 
; (-1) ; 
X46 O = lt 
; (=1) 05) on 
; (-1) (0) a 
(-1)(1) _ 
ola rows: X15 


This was the “key row"s 
New row: Ko 


old value —1 


BOW Ve US. ceva nanes) Go ein 


2) a cs 


Beet ©. sbter che new Cebleau is complercea aispece 
row Be ~ C5) to determine if there are any negative values 
les, Z,- c,< O); if all are positive or zero 
Ze = = < O) an optimal solution has been reached and no 
further iterations are necessary. If any Le “ oe < O an 
Optimum solution has not been reached and further iterations 


are required: one begins again at step 1, going through the 


entire procedure (steps 1 through 6). If, however, some 





G2 
e - S. moana. a iil a5 < 0 then the solution is unbounded 
or infinite (further iterations would be of no benefit since 
an optimum solution does not exist). 

Once an optimum has been reached the solution is 
found for the vectors in the basis (x, ) in their respective 
meng Under column (b). "Shadow prices" of variables not in 
the basis may be found under the respective columns (X 5) 
in row — - Cy). Shadow prices represent the solution of 
the dual problem. As mentioned earlier, it should be noted 
that differences exist among the various authorities rela- 
tive to certain aspects of handling the simplex computations. 
Paramount among the differences is that some authors do not 
negate the objective function initially when performing a 
minimization operation. This results in a change of 
criteria for selecting the variables to enter the basis at 
each iteration as well as a change of criteria for deter- 
mining an optimum final solution; for example, Greenwald 
follows this method. >” Consistency is an all important con- 
cern in using the simplex method, and these variations are 
of no great difficulty if one is consistent; however, one 


should be aware of procedural differences existing in the 


Ba cCerature . 





"Greenwald. Op, (Cie. PP - OS=607 
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The Modified Distribution (MODI) Method 


The modified distribution method is a somewhat less 
involved method of solution than is the simplex method. 
The MODI method, being applicable to only certain types of 
problems (generally called “transportation type" problems), 
is not a universal method of solution as is the simplex 
Somputcational procedure. 

The transportation type problems to which the MODI 


method may be applied take the general form, as follows: 


™m n 
Objective Function: s y c,. X,, = 2 (Min) 
a so 
Constraints: i 
55 . Dee a. JL eee eee 
jes fe 
al 
ye X,. = b De tas © 0 OD 
seh hee ae 
m n 
iS a. = » b 
i=l * jel 2 


HV 


Ks OF Wai 


These problems may be identified in actual situations 
Whewmweamparticiular system is in balance (i,e., when availe 
able resources exactly equal the requirements or demands 
for the resources), and the resource to be allocated among 
the various demands is homogeneous. In addition, the 


resources and demands must be expressible in a single unit 
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(such as gallons, pounds, hours, etc.). Realizing that we 
are still discussing linear programming with its inherent 
linearity requirement, it might appear that these additional 
limitations would severly restrict the usefulness of the 
MODI method; actually this is not the case. Transporta- 
tion type problems occur very frequently and the MODI method 
is often applied with great effectiveness. 

As before, we shall use an example to demonstrate 
how the transportation type problem may be formulated and 


how the MODI method may be used for its solution. 


(Transportation Problem) 


Assume that a U. S. Navy Supply Officer has requisi- 


—ELOns ELrom five destroyers, each ordering various quantities 


of a particular lubricating oil. These demands are as 
follows: 

Doo UR© Vk 1 2 3 4 2. 

DEMAND 

QUANTITY 100 60 40 jie. 29 

(GALLONS ) 


The supply officer has three storage locations where he keeps 
the lubricating oil and he knows that to get the oil to the 
various piers where the ships are tied up it will cost the 
government money as shown in the following transportation 


cost table: 
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COST (¢/gal.) OF TRANSPORTING LUBRICATING OIL FROM 
STORAGE LOCATION TO DESTINATION (DESTROYER) 


DESTROYER _1_ a ue Lou 5 

STORAGE 

LOCATION 
1 3 2 3 4 1 
2 4 1 2 4 2 
3 1 O 5 3 2 


The supply officer Knows that he will have no trouble fil- 
ling the requisitions because he has exactly sufficient 


BMierieating @il On hand to meet all of the demands, as 


follows: 
STORAGE OU AN TIT yeaGaa les) 
LOCATION AVAILABLE 
i 100 
Z IgA 
S| 7 


mie peoolem facing the supply officer, then, is to get the 
desired quantities of lubricating oil to each destroyer at 
the least cost to the government. 

The tableau format is the generally accepted compu- 
tational data display device used in the MODI method as it 
was in the simplex method; however, the MODI tableau is not 


the same as the simplex tableau. 
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By examining the first tableau (3.2a) it can be seen 
~enat: 

a. The quantities of lubricating oil available at 
the storage locations are entered under column (S); 

b. The requirements of the various destroyers are 
entered in row (D); 

@.  Transpertcation costs associated with shipping 
from a storage location to a particular destroyer are 
entered in the upper left hand corner of each cell. 

Gee Souantet1es Of lubricating Gil. to snip £rom ome 
storage location to a particular destroyer are entered 
within circles in the center of certain cells. The deter- 
mination of these quantities is explained below; in the 
first tableau the circled quantities represent a first 
feasible solution; 

e. Column (u) and row (v), which are left blank in 
the first tableau, will be utilized in later computations; 
these may be considered, simply, as two arbitrary variables. 
Results of computations with these variables will be 
emeerecd@in the lower richt hand corner of later tableaux. 

To begin iterations one must, as in tne simplex 
method, find a first feasible solution (i.e., a solution 


which is not necessarily optimal but is a possible solution). 





yah 

There are several ways available for obtaining a 
first feasible solution in transportation type problems; 
Hadley describes five different methods of getting the 
first feasible solution, ~° 

One effective way of getting the first feasible 
solution is to scan the tableau and find the smallest 
shipping cost (if there are "ties" where two or more costs 
are equally tne smallest, then anyone of the costs may be 
chosen arbitrarily). After the smallest cost has been 
determined one assigns as much as possible of the item 
Being transported to that cell having the smallest cost. 
One then finds the next smallest cost and again assigns 
as much as possible of the item being transported to that 
cell having the second smallest cost. This procedure is 
mOllowed (with the third smallest cost, etc.) until all 
destinations' demands are satistied and all sources are 
exhausted. For example. in tableau (3.2a) the smallest 
cost is zero (0) in cell (3,2) and we assign as much as 
Beet Dleteor the tobricating oll to. chis box (or 60 gal.) 
The next smallest cost is l¢/gal. in cells (3,1), (2,2) 
BMGs lle o)- we Gan assign 15 cals-(which is the quantity 
that can still be furnished from storage location 3 
L75 galt - 60 gal.]) to cell (3,1) and 25 gal. to cell 
(1,5), but we can assign nothing to cell (2,2) since the 


destination requirement (60 gal.) has already been filled 





40 
Hacleyes Oo¢. Cit. , s.°293) and pps 304=309;; 
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ae cell (3,2). The next (or third smallest cost is 2¢/gal. 
maecells (1,2), (2,3), (2,5) and (3,5). Immediately we can 
rule out boxes (1,2), (2,5) and (3,5) since the destination 
requirements have already been filled. This leaves only 
cell (2,3) to which we can assign a quantity of 40 gal. We 
continue in this manner and finally obtain the first feas- 
ible solution as shown in tableau (3.2a). 

We are now in a position to begin iterations toward 
a final optimal solution, in a step-by-step procedure. 

Step 1. The vacant cells must be evaluated to 
determine if the feasible solution is optimal. To do this 
we use the two arbitrary variables (u) and (v), making the 


following definition: 


B 
Co eee oa ey 
1) al J 
where : oa = the cost associated with the cells 


(i,j) in the feasible solution which 


have been assigned quantities. 


miei, HOrawany particular cell (i,j): 


Hadley gives a concise proof of this relationship. “> By 


arbitrarily assigning one of the variables (us, ver some 





41 Hadley. Op. Her eey, De. SLO. 








yes 


value we can solve for all Zz. For example, in 


oe Ge, ae 
1) 1) 


tableau (3.2b) we have arbitrarily set u, = 0, then: 


B —_— e -—- t — e — 
C14 = 3% Ws G-. then, 2: =-0 = Mae Vv, = 3 
oP Sean. == 38 Chen. ea =—U + 3° Uu = | 

Bal : il : : 2 : 2 
Ae a Sear =e Bere say a ep 

Sik ’ al : ‘ 3 : 3 

Gee ea OO en a ee 
So) ile ila | ae ! = Dy aoe 
oB =e = 1) ws le Ener: SSA Ce we en er 

23 : 2 : : Cy 3 
oB =e Are 1 ))2. = olive ee @ (=o ee 4 ae — ee Ag ee 

Baga : 2 ‘ : 7 AN 4 

ae Me A ie Oe VS Wedges 1 ete ee ye ee PL 
ea ca : aL Es r 7 = 5’ 5 r= 


Once all us, and see have been determined, as shown above, 
we may evaluate the empty cells. For example: 


Z, Ou se2 Ge a7 eH 


+- 
< 
| 
Q 
I 


in Oe 


| 
< 
a 
< 
| 
0) 
I] 


z Ss 


33 (-2) + 1-5 = -6 


e30n 
This procedure is followed until all empty cells are evalu- 
ated: the resulting evaluations are entered in the lower 
right hand corner of the empty cells in the tableaux. 

pies Pec a Ener tableau as next scanned tO, determine ar 


any Z. 


ij 7 oan is positive. When any Z, PO San 


ii ae aks 


A 


Optimum solution has not been reached. If all ae - Cij 


an optimum solution to the problem has been achieved. It 
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may be noted here that this criteria is Just the opposite 
of that used in the simplex method. Some authors use a 
somewhat different procedure with the resulting criteria 
being if any evaluation of an empty cell is negative then 
an optimum has not been reached: (for example, see 


cae As long as one is consistent in the evaluations, 


Garvin 
the criteria is unimportant. The difference is pointed 
out here only as a caution that there are also variations 
Se this procedure among the authorities, as there were 
differences in the simplex method. 

Step 3. When an optimum solution has not been 


reached the empty cell containing the largest positive 


value Of Z..-=- cCc.. 
1) LJ 


is chosen to be filled by some quantity 
presently filling some other cell (if "ties" occur among 
the ereeee value of ae - os any of the largest positive 
values may be chosen arbitrarily). Since the system must 
be in balance at all times, several adjustments must be 
made. The most efficient method of making these adjust- 
ments is as follows: 

a, Place some identification (say: theta [9]) in 
Ene cell to be filled (notice cell [2.2] in tableau [3:s2b]}). 


b. If we add some amount 9 to a row, then, to 


remain in balance the same amount 9 must be subtracted from 





42carvin. Oe Gd ete Dey soe 








Tie, 
Bome other circled quantity in that row: but, from what— 
ever column in the row we subtract 90, then, we must add Q 
to another row in that column; etc. In other words, we 
Begin at the cell having the largest positive value of 
Sag ~ Cis and add an amount 9 to that cell, then, we pro- 
@eess from that cell alternatively subtracting and adding 
9 in cells having circled quantities until we return to 
the cell at which we began. For example, in tableau (3.2b) 
we place a (+0) in cell (2,2); a (-8) in cell (2,1); a 
(+0) in cell (3,1)3: and a (-90) in cell (3,2). Thus, we 
have completed the circuit. Notice that this route which 
has been followed is the only one possible. If we had put 
a (-9) in cell (2,3) instead of cell (2,1) there would have 
been no other cell in column (3) to which we could have 
added a (+0). 

c. Once we have determined which cells have to be 
adjusted we find some cell having a negative adjustment (-@) 
and which has the smallest circled quantity among all the 
negative (-@) adjustment (in tableau [3.2b] it is cell 
Pel eee Eheta (0) is then set equal» to this: emallest cirelied 
value and the respective additions and subtractions are 
made, with the results entered in corresponding cells of a 
new tableau. Cells which were not affected by theta (0) 
adjustments are reproduced in the new tableau from the 


previous tableau in their respective positions. 
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One iteration has now been completed and to continue 
on to reach a final optimal solution we must return to 
step 1, and repeat the entire procedure again until all 


t 


vacant cells have Z..-c.. 
bey) ee 


= O. In this example problem 


we reach the optimal solution after two iterations (see 
maopleaux {[3.2c} and [3.2d]). 

To interpret the final results we merely examine 
the final optimal tableau and read off the various quan- 
tities to be shipped from a source to a destination. In 


the example problem we see that the source S. (storage 


Z 


Mecation 2, in row 2) ships 60 gallons of lubrication oil 


Eo destination D. (destroyer 2, in column 2). It is inter- 


esting to note that storage location S. gives nothing to 


3 


destroyer D. even though the vessel must be tied up at 


Z 


S3's pier because there is no transportation cost involved 


211 Shipping from S3 fe. D.- 

To obtain the cost of a shipment one merely multiplies 
the cost ic; 3) associated with a cell and the quantity to 
be transported (circled quantity in the cell). The total 


cost of the optimal program is, obviously, the sum of all 


Ene individual shipping costs. 
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TABLEAU 3.2a 





Column >l 2 3 4 5 


g Kon 
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TABLEAU 3.2b 


~ 
in 


Fas 


(Fo ese) 
(22 Cc) 
(10)(4¢) 
(40)(2¢) 
(75)(4¢) 
(1s )Cle) 
(60) (O¢) 


Total 


G75) 3c) 
(25) <1¢) 
(10)(4¢) 
(40) (2¢) 
(75)(4¢) 
fee Jee) 
(60)(0¢) 


Total 


TOTAIS COST 


t G 
N be b- 


TOTAL COsSe 
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$2.25 


OS. 
0.40 
0280 
i 018) 
Oz 


0.00 
$6.85 


$2.25 


OR 25 
0.40 
O250 
3.00 
O25 


0.00 
$6.85 



























aL: 2 a 4 oS S., u 
eo umn > a a 
ROW —1 2 3 
* as) 
O O 
4 at 2 
5 
3 
-4 
TOT: COST 
De 40 (75) (36) = 
| Cs (Ok ae — 
(160) (17) = 
(40)(2¢) = 
‘5 3 (75)(4¢) = 
(25) ie) == 
(50) (0G): = 
Total 
TABLEAU 3.2¢c 
Column pal Z 3 4 D Ss; ee 
—_ 
l (25) 100 O 
1 
25 -l 
O 
2 
73 -2 
O -l 





BOOS 


| (50)(4¢) = 
O75) 2 Gilet gee 
3 (60)( ie) os 
° : CAO 2) c= 
PGs) — 
(7S) 4a) c= 
Total 


Final (Optimal) Tableau 


TABLEAU 3.2d 


LOTTA. Cost 
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$2.25 
0.25 
0.10 
0.80 
3.00 
0.25 


0:00 


$6.65 


SO- 7S 
2500 
Os 25> 
0.60 
0.80 
i200 


On 7D 


$6.15 
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The "Stepping-Stone" (Transportation) Method 


What is generally called the transportation (or dis- 
tribution) method is associated with the "stepping-stone" 
technique (which is another method of evaluating the vacant 
cells in the transportation tableaux). This one difference 
distinguishes the transportation method from the MODI 
method. Most authors of linear programming texts give 
comprehensive treatments of the “stepping-stone” method 
(for example, see Charnes and pooner 

The transportation method, using the "stepping-stone" 
technique and the MODI method, using the evaluation techni- 
que described above, obtain the same numerical results. It 
is believed, however, that the "stepping-stone" method is 
more difficult to work with and is more apt to cause a 


person to make more calculating errors. 


Degeneracy 


Degeneracy may be encountered in either the simplex 
algorithm or in thé various distribution algorithms. 

In the simplex method, degeneracy is caused by more 
than one vector passing through a single point, and it is 
manifested when, in trying to determine which variable to 
eliminate from the basis, “ties" occur among the smallest 
ae 


43 
Chaknes ana Cooper. -Op. Cites, iep. 41-49. 
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e 

rwratios. In the case of “ties", regardless of which variable 
is selected to leave the basis, the other variables involved 
in the tie will be reduced to zero, and there will be no 
change in the functional (Z). In general, the procedure 
for breaking ties, as given above, will be sufficient to 
allow the procedure to continue to a final optimal solution. 
It is possible, however, that the computations will cycle, 
or, in other words, after several iterations, some pre- 
viously eliminated variable re-enters the basis. It should 
be obvious, then, that cycling and degeneracy are not 
synonymous. Degeneracy is necessary before cycling can 
occur, but degeneracy is not sufficient to cause it. Garvin 
states that ". . .while degeneracy is common, cycling is 
extremely rare."*4 Although there have been problems 
devised to cycle, Gass reports ". . .there have not been 

any practical problems that have been Known to cycle"? 
Since degeneracy and cycling in the simplex method are of 

no practical concern the method of solution for the simplex 


algorithm described above would appear sufficient; how- 


ever, there are procedures designed to absolutely prevent 





“4garvin. Op. Cit., p. 47. 


45 
Gacsy saul, a. Linear Programming Methods and 


Applications, New York: McGraw-Hill Book Co., Inc., 1958. 
Bee OSs 
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cycling and these may be found in the literature (for 
example, see: eae, Garvin, */, or padiey. 

Degeneracy may occur in the various "transportation" 
methods whenever there are less thanm +n -- 1 cells of the 
tableau filled with circled quantities (m = number of 
Pources Or rows: n = number of destinations or columns). 
This situation may occur at any stage of the computation. 
and may be overcome by a very simple procedures: one simply 
fills sufficient vacant cells with zeros in positions which 
will prevent degeneracy. The cells containing zeros are 
then treated as quantities to be shipped (circled 
quantities). Hadley presents a similar method which pre- 
vents degeneracy in transportation computations, from the 
eutset. 7 

It is theoretically possible for a transportation 


type computation to cycle: however, it would appear that 


cycling is even less of a problem, here, than in the sim- 








plex algorithm. Hadley states: ". . .no transportation 
problem has ever been known to Clones 
Sorbian. Ghapter 7. 
47 ; 
Camvilic, Op. Cll. , Chapter, 14. 
48 


Hadley | Ops Cie.) Chapter 6. 


emesis. ton. SOR 


erga eie HCN. 





CHAPTER IV 
GENERAL APPLICATIONS OF LINEAR PROGRAMMING 
General 


Examples of linear programming applied to various 
problems of industrial management are legion. Perhaps one 
of the reasons for such wide-spread application of eects 
programming in industrial enterprises is the relative ease 
@emdatce Quantification such as costs, profits, mix ratios, 
eee. Im addition it is a generally recognized fact that 
for management decision making, complete data accuracy is 
not always imperative; estimation of data is, in many 
cases, acceptable and adds, therefore, to the relative ease 
Sr tne data quantification. 

In this chapter a survey of the generalized classi- 
cal applications of linear programming will be presented 
together with five specific numerical examples. It is 
intended only that this should indicate the manner in which 
linear programming may be applied, and it is not intended 
to be a complete representation of all the applications to 


which linear programming may be put. 


Classica l Applications of Linear Programming 


There has been a tendency in the literature to 


Classify applications of linear programming by means of 
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sterotypes. In one sense there is nothing wrong with such 
meorocedure;: in fact it aids classification of applications 
considerably. On the other hand, however, when an idea, 
concept or application (as being discussed, here) is stero- 
typed, it is generally very difficult to remove the idea, 
concept, or application from one viewpoint and put in into 
a different perspective and setting. The real worth of 
ideas, concepts, and applications, generally, lies in their 
versatility. 

In this section some of the sterotyped classifica- 
tions of linear programming applications will be pointed 
out, together with some possible diverse uses of the par- 


ticular applications. 


(Diet Problem) 


The first classical linear programming problem deals 
with the feeding of a group of individuals (say, Navy 
personnel aboard ship) where it is desired that the per- 
sonnel be given meals meeting all nutritional requirements 


but with food of minimum cost. Thus, if we let: 


Xs = 1 th food item in the diet 

c; = cost of food item x. 

aay =) Mubricional content. (OL nutrient b,) ala 
food item xX. 

b., = j th nutrient required ina diet. 





A nen : 


Se jective hunecien: 


Constraints. 


Although a linear programming 


ning will indicate how to feeda 
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et =e LOM aa) 
> ‘ 
ifi = 2; et areal re * el, 
“Si 20 (i=1,2,3,. ..,n) 


approach to meal plan- 


group of Navy personnel 


with food containing the required nutritional elements at 


minimum cost to the government, 


it is not a generally 


recommended method of general mess management since it may 


lead to some rather unpalatable meals. 


being very objective and single-minded, 


Linear programming, 


Gan. InGil@ave Mot ming 


about food variety, nor can personal food preferences (among 


Navy personnel) be optimized. 


The diet problem was the forerunner of all the blend- 


ing or mixing type probiems of which the widest application 


has been in the petroleum refining industry relative to the 


Blending of gasoline. 
blending problem is given in the 
chapter. 

The diet problem may also 
such as the manufacturing of any 
sausage to paint or steel) where 


CUCtCION s1G Gesirea and there are 


A numerical example of a gasoline 


rellowing §<section of hve 


be applied in situations 
product (from ice cream or 
a minimum cost of pro- 


certain minimum raw 
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material requirements: or, restrictions on the availability 
of raw materials: or, restrictions on the quality of the 
Peeauct. Oh, FreSurIecei1ons On Che quantities: of CUtpUE OF 
mhe product which may be produced; etc. Generally these 
problems take the form of having an infinite number of ways 
in which the raw materials can be combined to make the out- 
muc product, 

Blending, or mixing, type problems are not limited 
to inanimate product manufacturing. Another aspect of 
blending was described in chapter II when the output of a 
disbursing clerk was maximized and his time was "blended" 
or “mixed” in working on pay records and vouchers. Simil- 
arly the utilization of equipment could be handled in this 


way. 


(Transportation Problem) 


The transportation problem was the second type of 
linear programming problem to be specified and, perhaps, 
is one of the most, if not the most important linear pro- 
gramming problem type. Its importance lies in the fact 
that a relatively simplified algorithm is available (the 
MODI method) for solution of transportation type problems 
and has a great range of applications. Since the general- 
ized form of the transportation problem was given in 


chapter III it will not be repeated here. 





86 
One application of the transportation problem was 
indicated in chapter III; that of simply finding the least 
cost of transporting a homogeneous product from several 
storage locations to several destinations where the total 
demands for the item, at the destinations, exactly equalled 
the total availability of the item at all the storage loca- 


tions. There are several variations of this one applica-. 


iON $ 

a. Availability of the item exceed the demand for 
it: 

b. Demands for the item exceed the availability 
eo. it, 


c. Only partial quantities can be supplied from one 
storage location to some destination. (A numerical example 
®t this variation is given in the following section of this 
enapter. ) 

d. Items can, or must, be transhipped between some 
destinations. 

Another application of the transportation problem, 
of importance in the Navy, is the routing of sea going 
vessels (say fleet tankers). In this application the 
tankers are located at various ports, anywhere throughout 
the world, and they are required at some other ports to 
transport fuel to some other locations. The problem, 
obviously, is to find the optimum assignment of tankers to 
load the fuel and transport it at minimum cost to the 


government. 





oy 
Ojewmoruher application o= —Eme eransportation problem 
is its use in bid evaluation. Here, several bidders pro- 
pose to supply the government (at various government storage 
locations) some required item, at various costs. By utiliz- 
ing the transportation type problem approach, assignment 
of contracts to the various bidders may be made at least 


cost to the government. 


(Caterer Problem) 


The classical presentation of the caterer problem 
concerns the manager of a restaurant (for example), or a 
caterer, who knows that on specific days in the future, 
while serving meals, he will have various requirements for 
napkins. To begin his business he will have to buy an 
initial quantity of napkins but thereafter, the manager, 

Or caterer, may either buy more napkins or send the soiled 
ones to laundries which can return clean napkins in speci- 
fied times at specified costs (with faster service naturally 
costing more money). The problem, then, is to find the 
Optimum program which will give an adequate supply of nap- 
kins at all times at the least cost. 

The generalization of this problem may be accomplished 


as brolvows: 
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Let: 

R , = elean napkins required on the j th day ( j=1,2,3. .,n) 
p = days required for slow laundry service 

q = days required for fast laundry service 

2“ P 

a = cost of new napkins 

b = cost of laundrying one napkin with slow service 

cS = cost of laundrying one napkin with fast service 


Ess Cc > b 


x; = napkins purchased on the j th day 
y. = napkins sent to the laundry with slowest service, on 
j 
the j th day 
Z. = napkins sent to the laundry with fastest service, on 


J the j th day 


S = soiled napkins stored on the j th day and not sent to 
any slaumary . 


Then: 


and 


The objective function is then: 


n 
5 (ax. + by. + cz.) = Z (Min 
a j ve 5? ( 


A numerical example (presented as the Wardroom Steward 


problem) is given in the following section. 
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The caterer problem may be applied to procurement 
and maintenance problems. For example, if a supply officer, 
in order to meet future specific demands, has to determine 
Pmether a certain type of spare part should be either pur- 
chased or repairable items sent to various repair snops 
(assuming other repair criteria has been met) then we have 
an analogy between napkins and spare parts. 

The storage parameter (s 5) may be assigned some 
cost, also, and then the caterer problem may be used to 
determine the program to optimize over-production (or, 
excess procurement) in one period (with attendant storage 
requirements) and future demands for some item. In this 
manner production (or procurements) may be "smoothed" when 
there are seasonal fluctuations in requirements for the 


item, 


soument On the Comouter Solution of Linear Programming 


Problems 


The numerical examples contained in the following 
sections were solved by means of an IBM 16020 computer, at 
the University of Kansas. The program used in obtaining 
solutions was from the 1620 General Program Library, titled 
Linear Programming II (10.1.008), by F. W. Wood. The pro- 
gram is a five (5) phase program written in "Fortran with 


Format" language. 
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Although there are published instructions for the 
program which was used (10.1.008), these instructions are 
dependent upon the published instructions for the program 
titled Linear Programming I (10.1.007); both sets of 
instructions are necessary to successfully work with 
Linear Programming II (10.1.008). 

The program negates the functional and, therefore,. 
actually solves a minimizing problem. To solve a maximiz-— 
mac problem the objective function first has to be 
negated, prior to solving. 

Since the library program instructions are not 
specific regarding the times required for solution of 
problems it may be of interest to note the times which were 


required for solution of the numerical examples which 


EOollow: 
Approxi- 
No. of No. of Itera- mate Time 
Vari- Equa- eLOnsS Required 
Problem ables EROS. eRe CUtme cl CMG.) 
Gasoline Blending ZZ 10 14 10 


Pa stribution (Bal= 
anced supply & 
demand ) 40 ea 36 fs 


Distribution (Unbal- 
anced supply & 
demand) 44 LS 40 23 


Wardroom Steward 40 20 2S es 


Wardroom Steward 
(Moetaieation. now 1) 50 Zo 26 18 
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Numerical Examples of Linear Programming Problems 


This section contains five numerical examples from 
the general classifications of linear programming problems 
presented in a previous section of this chapter. Although 
the numerical examples are somewhat more involved than have 
been presented thus far, they are considerably simpler than 
they would otherwise be in real life situations. eyo 
theless, the following examples are intended to indicate: 
(a) how a problem may be put into the farmework of a 
mathematical model; (b) how the mathematical model may be 
mee into a proper format for computer solution; and 


(c) finally, how the computer solution may be interpreted. 
Sak 
(Gasoline Blending Problem) 


A company wishes to make as much income as it can 
from its refinery operations. [It is believed that the manu- 
Facture of aviation gasoline will be profitable. It is 
proposed that three types of aviation gasoline (types A, 

B, and ©) should be made, In order to do this, raw 
materials consisting of: alkylate, catalytic cracked gaso- 
line, straight-run gasoline, and isopentane, are to be 


blended to make the aviation gasolines A, B, and C; any 





eta ea@apeed vem Gale. FO. Le pow Of-Co- 











ge 
excess of the materials will go into gasoline for automobile 
use. The properties of the above raw material components 


are given below in table (4.1). 


TABLE 4.1 


GASOLINE BLENDING PROBLEM 
RAW MATERIAL COMPONENTS 


OE ON TE 


Octane No, Soto eNS 
Reid | rola s| : - | Availa- 
















(with oe 
Vapor | 0.5 ml/gal. biliey 
Oe e Pres- | Tetraethyl-1{2 ea (bbls/ 
Tetraethyl- 
sure lead con- day) 
lead) 
tent) Pare | 
Alykylate 107, p07. Ses 
Catalytic cracked | Seraw 27,052 
gasoline 
Straight run gaso- 4 5A 
Pas 





12 
ee oa 


Lsopentane ; 2025 eaten els 


The raw material requirements and expected revenue 





Of ene Various products are given below: 





a3 
TABGE, 42.2 


GASOLINE BLENDING PROBLEM 
RAW MATERIAL REQUIREMENTS 


Reid 
Tetraethyl Income 
ressure 


LOO 6.042 


Av. Gas. B 


arcane | £2 
awtoces, | | =e < 


The blending process is shown in figure (4.1). 





lV 


The following constraining relationships may be formed by 
mBeilizing the data found in tables (4,1) and (4.2): 


Reid vapor pressure constraint: 


AV. 5X1 + 8X. + 4X. + 20.5X, = T(x) + X, + xX, + Xy) 
GAS. z= 
A Or 2X1 + X5 = 3X3 + 13.5X, —ae 8) 
— 
AV. 5X, ex. t 4X, + 20.5X, = 7(X, + Xo + X, + Xg) 
GAS. = 
B Or ~2X_ + Xe - 3X5 + 13.5X, =" 
< 
AV. 5Xg + 8X,9+ 4X,,+ 20.5X,5= 7(Xo + Xyot X14+ X15) 
GAS, Z 
C GE -2Xo + X10 — 3X, ,+ 13.5X)5= 0 
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ALKYLATE i . 
ScArTALYT LC 
CRACKED xX 
PASC OT ; Z : 
“ AVIATION GASOLINE 

STRAIGHT RUN 

D4 AN it 
GASOLINE eg) = 


ISOPEN TANE LL xX, 


Vey 





AVIATION GASOLINE 
a Sa 







AVIATION GASOLINE 
he ke 





Wy OW Wy 


AUTOMOBILE .GASOLINE 





FIGURE 4.1 


SCHEMATIC DIAGRAM OF THE GASOLINE 
BLENDING PROBLEM 
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Octane Number constraints: 


> 
AV. 94x, + 83X, + 74x, + 95x, = 80(X)+X,+K +X, ) 
GAS. . 
A or 14x, + 3X5 — OX, + 15x, =O 
> 
AV. 107.5X,+93K.+87xX_ + 108X,, = 91(X,+X_ +X +X, ) 
GAS @ > 
"BR or 16.5X,+ 2X67 4X + 17X, =) 
> 
AV. 107.5X9+93X, 9+87X, , t108xX, , = 100(X9+X 9 tX1 4 +X) 9) 
GAS. = 
tt 40 oes — — 
Cc ONG 7.5Xy 7X16 13X,,+ 8X, 5 O 


Raw material availability constraints: 


IA 


at eo ta 


1 5 9 7300 


IA 


X~ + X. + X 


2 6 Noes cee 


A 


> ee aa? Ome on 


fe 7 ea 


IA 


Doet? eat OK 


4 g (eeeo ee 
The total volume of liquid is 11,833 bbl/day (3,800 + 
3,652 + 4,081 + 1,300). If we let the respective volumes of 


aviation gasoline A, B, and C which is produced bes: 


Vai Var Vas 

Then: Va = xX) ~ X5 + X53 + Xy 
Va = Ke oa Xe + Xx + Xp 
Von = he or Xa et OX 


c 2) 10 te 12 





36 


The automobile gasoline produced (Vaan) is, thens 
a 


eae OS 3 (Vn 


AG fever 


B 
The objective function to be maximized, then, is: 


4.908V,+5.437V, + 6.042V,,+4.548(11,833-V,-V,.-V,) = Z% (Max.) 


Or ¢ 


53,816,484 + 0.360V, + 0.889V., aa 1.494V, = Z (Max. ) 


Temporarily the constant (53,816.484) may be ignored and we 
have the objective function and mathematical model as 


follows (after insertion of appropriate slack variables): 


OBJECTIVE FUNCTION ; 


0.360(X, + Xo + X4 + X,) + 0.889(X, + Xe + Xo + Xa) 


ne 1.494(Xg+K, 4 tXy +X, = Z' (Max, ) 


>) 
Note that the objective function shown above is negated when 


preparing the input data for the computer program. 





CONSTRAINING EQUATIONS : 


_ 2X, + Xo _ 3X3 + 
14x, + 3X. - 6X. + 
_ 2Xe + Xe _ 3X2 + 
16.5X,+ 2X— ~ 4x - 
~ 2Xo + X10 - 3X1 ,+ 
7.9%9- 7X10 - 13X,,+ 


13.5%, 
is OX 
1S ie pA 
17.0X, 


13.5X)5 


Ls 


14 


15 


iG 


tUY 


ic 


oy 


II 
O 


= Q 


= 3,800 


= 2,652 


= 4,081 


= 1,300 


The constant which has been ignored in the objective func- 


tion represents the income which would be 


kK. = O? Or, in other words if no aviation 


J 


made at all. After solution of the above 


be remembered that the constant has to be 


Picome indicated for the functional since 


Obtained if all 
gasoline were 
system it must 
added to the 


the functional 


Will indicate only the additional income derived solely 


from the manufacture of the aviation gasolines. 


Diagrammatically the solution is shown below in 


migure (4,2). 
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ALKYLATE %, = 0-00002222 

CATALYIACr = 

CRACKED | _ 

ee ee ee | SAVIATION GASOLINE 
STRAIGHT NAM 

RWNonine |[%3] = 9-90095185 | (0.00007408 bbis/day) 


Ox. 
ISOPEN TAN Led - “i | 








34. 7 


AVIATION GASOLINE 


a 
i 2961. 110 (3658.359 bbis/day) 


6632 50m 
S/ Ga 40 Oke. 









a | 
uy 2652.0001..' AVIATION GASOLINE 
7 Uy Oe 
1119.8894 
il) 8 ! sal (8147.5474 bbls/day) 
NM 610% 257.0.) 
LA i) eee 


AUTOMOBILE GASOLINE 


(26.6 bbls/day) 


PLGURE, “42 
SOLUTION OF THE GASOLINE BLENDING PROBLEM 





w 


eh 

The functional, as given in the computer output, is 
$15,425.161; adding the constant term $53,816.484, the 
total income, then, is $09,241.645. It can be noted, how- 
ever, that the solution calls for a minute quantity of 
aviation gasoline "A", to be produced, which would be 
neglected in reality, as being an unprofitable product. 
manegiecting all production o£ aviation gasoline “A" our. 
functional remains the same, however, since the small 
quantities called for in the optimal solution add a 


$0.000024 and cannot be detected in the final solution. 


(Distro tien Problem)? 


(Balanced Supply and Demand) 


AU. S. Navy Supply Officer has an item (for example, 
say cement, but the product is really immaterial) which is 
stored at four depots. He receives requisitions (or 
demands) from ten regquisitioners (or customers) for various 
quantities of the cement and the officer Knows what shipping 
costs are involved when shipping cement from any of the four 
depots to any of the ten requsitioners. These data are 
shown in table (4.3). Each cell of table (4.3) must be 
described by a variable because each cell is associated 


with a specific shipping cost and represents some potential 





> Data adapted from: Hetrick, James C., "“Mathemati- 
cal Models in Capital Budgeting," Harvard Business Review, 
(January~February, 1961) reprinted in HBR Statistical 
Peciston Series, part IL, p. 67- 
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suiantity of cement to be shipped from a depot to a requisi- 
tioner. In the final solution the variables will indicate 
which depot should supply which requisitioners and in what 
quantity. Then, variables may be assigned as shown in 
Gable (4.4). 

Recalling the generalized transportation problem 
formulation shown in chapter III, a mathematical model is 


Peonstructed as follows: 


TABLE 4.3 


SHLEPEING COSTS OF “THE DISTRIBUTION (PROBEEM 
(BALANCED SUPPLY AND DEMAND) 


eae Storage 
Shipping Costs ($/cwt) vail.| Location 
(cwt) | (Depot) 





Meet 1 .15)0.21)/2.5411.98)/0.49) 1.001 0.36 





262 oe OUlyoD00 if 


1.67]2.45] 0.50] 2.80 2230 20 1.26|0.75 2a So S25) 32 Oe 


eS} 1.20| 0. 30) Zee OOO eGo. OTD 55 2.40 }3.00} 4800 


o.25] 0.65] Oa 2% 00|.25|0. Sonu. we a oe 15 [2.25 4800 a 
Ova nt ley 
2000] 2400! 1500/1300 | 1800;1700) 2200/1800} Requisitioned 
Cw) 











TABLE 4.4 


VARIABLES OF THE DISTRIBUTION 
(BALANCED SUPPLY AND DEMAND) 


VARIABLES 
x, X5 K 3 Ky Ae LX 
Ko Kos 





“on 


























ees 


OBJECTIVE FUNCTION : 


ak, + 1, 15X.) + O.21X 


a 


2 


0.36X, Pec eooKe + 13, 00X 


2) 


2.15X,-+ 0.80X 


15 ihe: 


mh. 20X,,+ 0.30XK..+ 2.99X.~,+ 


Ze Zo 


2.40X. 4+ 3.00X 4+ 722% 


0.38X,,+ 0.89X,.+ 0.25X,4+ 


50 oe) 


3 


10 


+ 1,260X, 5+ 


SL 


+ 3.54X 


+ 1.6/X 


Oo om 


2.00X 


+1) 6px 


Le fo 






1500) 300 53001 2100) £800] 17001-2200 


4 


gle 


igs: 


Z> 


32 


Sie, 


7 OS 


+ 1.45% 


+ 2,/5X 


+ 0,60X 


+ O.,10X 


+ 2eZzon 


PROBLEM 


at 


Xo 


5 


A 


19 


26 


a5 


40 


0 


oe 


+ 


ot. 


te 


oe 


EO 


Qty. Storage 
Avail.}| Location 
lew) Depot ) 


6500 


2000 


0.49X- we 


0.50X, 3 


3.25X56 


1,10xX.- 


2.00X., 


Zhan) 


+ 2.80X 


telex 


+ 125% 


1 


Z 


Quantity 
s{ 7] ° os 


1,O0O0X. + 


147 


Dae 


+ 0.99X5 4+ 


si 





GZ 
Note that this objective function does not have to 
be negated prior to preparing the computer input data; since 
the computer program is written to negate the objective 
function automatically it, then, causes the computer to 


eOolve for the minimum. 


Constraints: 


Supply (or row) restrictions: 


Pe eee 0,200) 


is 
va 
I 


2000 


M 
ms 
I 


4,800 


xm xX. = 4,800 
j=31 J 





1HO)S) 


Demand (or column) restrictions: 


Xy ee X14 + X54 + X35 — ee O18) 6. 
Xo + X15 + X55 + X25 = 2,400 
X3 + X13 + X53 + X33 Zz be SUO 
xX) + Xia + X54 + Xa, = OO 
Ke + Xis + X56 + Kae = 1300 
Xe + X16 + X56 + X36 = 27,190 
Xo + X15 + X54 + X37 = 1,800 
Xe + X18 + X59 + X39 = 1,700 
Kg + X19 + X59 + X39 — awe) 8. 
ea ok os ar = 1,800 


10 20 30 40 


The above system is shown in the canonical form of 
linear programming statement and is suitable for submission 
to the simplex algorithm rather than the transportation 
eeigorithm, 

The solution for this problem is shown in table (4.5), 


below: 
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TABLE 4.5 


SOLUTION OF THE DISTRIBUTION PROBLEM 
(BALANCED SUPPLY AND DEMAND) 


| 
| 
! 


Depot 


_— 
pel [see] [Tose Yr el Toe 

ne aff 
So LBOOT SOO £001 1 S00 epee] — 


ebebel Delo] es [ie [rewire 


© TOTAL COST: $23,113.00 








It is interesting to note that not all of the requisi- 
tioners are supplied by the least expensive method of trans- 
portation; for example, it seems that it would be less 
expensive to supply requisitioner No. 3 from depot No. 4, 
or customer No, 6 from depot No. 4. It is to be emphasized 
eaac the solution, as given above, optimizes the system, Not 
individual shipments. To optimize shipments to individual 
requisitioners, or shipments from individual depots would 
lead to sub-optimization of the whole system and, therefore, 


require greater overall cost. 





eS 
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(Unbalanced Supply and Demand) 


Suppose, now, that in the previous distribution 
problem the carrier who transports the cement from depot 
Mo. 1 to requisitioner No. 6 can handle only 1,000 cwt. 
because some of his equipment has broken down. (The 
previous solution calls for a delivery of 2,100 cwt. from 
depot No. 1 to requisitioner No. 6.) Since the requisi- 
tioner needs his material immediately, the question is how 
can the demand be filled, right now? 

To solve this problem we can proceed as before but 
mars time we must limit the quantity going from depot 
Mon lb tO requisitioner No. © to 1,000 cwt. To do this we 
change tre demand in column 6 from 2,100 cwt. and add one 
additional column having a demand of 1,100 cwt. The new 
column will have, in rows 1, 2, 3, and 4, the new variables 
X 


and X respectively. The new variable 


AT Si) G weer Kee 4a! 


Xa (a partial shipment from depot No. 1 to requisitioner 
No. 6) represents the quantity of cement which the carrier 
cannot handles therefore, in order to make sure variable 
Kay does not enter into the final solution it must be given 


a prohibitively high cost (say $10.00 per cwt.). Variables 
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x and X are identical, respectively, to variables 


Ao 13° a4 


X16: X56; and X36. therefore, they must have the same 


mransportation costs. 
Now the problem is the same as before. To formulate 
the model we simply add the new variables as part of the 


previous objective function: 


+ 0.38X aie Mia) 


Z + 10.00X fo. OO AA 


+ 0O.00X 


41 42 43 


Then, to the supply (or row) restrictions we must add the 


new variables, as appropriate: 


}4 


O 


- X, + X41) = 6,500 


J 


20 
y A. + X 
jeilan 


i 


2,000 


30 


ee eT OG 


a 43 = 4,800 


tf 


40 
a x. oh Oe =—a7 sod 
j=H3l1 
One additional demand (or column) restriction must be added 
and the restriction on column 6 must be changed from 2,100 


to 1,000, while all the other colum restrictions remain 


unchanged as before: 





LOG 


Xx + X + X + X =e OU 


Xai + XA5 + X43 + Xana =e Hae 


The solution of this new problem is shown below in table 


(4.6): 
TABLE 4.6 
SOLUTION VOF THE. DISTRICGLTION PROBLEM 
(UNBALANCED SUPPLY AND DEMAND) 
Depot 
feof — 06 
: 


2400] 1300 ie vs ae ae 1100] 4800} 3 
Le dr 
i700} | | _ 2300 = ) 1800 4800 
2000! 2400] 1500}130011300 1000 [1800 11700 | a se 
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Lt can be seen that this new problem has increased 
Beansportation costs by $24.00 ($23,137 = $23,113) but 
under the circumstances it is the optimum solution. 

It is interesting to note, in this problem, as 
before, not all of the requisitioners are receiving sup- 
plies by the least individual cost, and that by changing 
the route of supply to requisitioner No. 6, some of the 
other supply routes also changed. 

As mentioned previously this problem can be extended 
to cover other situations which occur in distribution 
eystems. For example, it 18 not imperative that the 
supply source availability exactly equal the demands from 
the requisitioners. 

“If availability of the resource at the depots exceed 
meouirements by the requisitioners then an extra column 1s 
added, as above. The extra column represents a fictitious 
(or dummy) requisitioner who is requisitioning all of the 
excess from the depots. The total for the dummy column, 
then, will be the difference between the total quantity 
of all depots and the total quantity being requisitioned 
by all the real (or actual) requisitioners. Since the 
dummy requisitioner does not exist and the excess material 
will not be moved from any depot the shipping costs assigned 


to the variables in the dummy column must all be zero. 
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Similiarly, if demands from requisitioners exceed 
availability of the resource at all the various depot, then 
an additional row (dummy depot) may be added with an 
“availability” equal to the total deficiency of ail the 
actual depots. Shipping costs associated with the dummy 
depot must be exorbitently high to prevent any "bogus" 
shipments from the dummy depot entering the solution until 
all the least expensive modes of transportation have 
exhausted the actual available supply. "Bogus" shipments 
will appear in the solution and will represent quantities 
of the resource which cannot be shipped to meet certain 
demands; for this reason the functional obtained in the 
eemputer solution will not be valid and an actual shipping 


cost must be computed by hand. 


(Wardroom Steward Problem)?" 


The Chief Steward aboard a large Navy vessel just 
being commissioned is faced with an optimizing problem. 
The wardroom has not yet been used and no napkins are on 
board. The Chief knows that the shipboard laundry will not 
be in operation for ten more days and that during this time 
meals will be served in the wardroom and he knows his nap- 


kin requirements will be as follows: 





> 4 ata adapeed. Pron Garvin. “Op. Cie. , pose Eo izo. 








HE DEG! 

DAY i 2 3 4 5 6 a 8 9 10 
NAPKIN 
me QOULREMENT 10) 60 30 70 318) 60 90 80 50 ieee) 

The Chief Steward knows he can purchase napkins for 
10¢ each (a = .10), and one laundry will launder soiled 
napkins for 1¢ each (b = .01) but it takes three days 
(p = 3) to get them back; another laundry will launder the 
soiled napkins for 3¢ each (c = .03) and return them in 
two days (q = 2). The chief steward must determine what 
mrogram to follow in meeting the napkin requirement at a 
minimum cost to the wardroom officers' mess, 

Recalling the problem generalized formulation pre- 
sented in an earlier section of this chapter, we can con- 


struct the mathematical model as follows: 





ces leah 


Napkin usage equations: 


Napkins Napkins 
Napkins Returned Re tu Gne a Napkins 
Pete UL Cchasea Giarees) eh ee (Slow) Required 

(x) ar (z) 5 (y) = (R) 
- x) = 50 
2 xX, = 60 
3 X 5 + X14 = 80 
4 Xx, + X15 = 70 
5 Xe + X13 + X55 = 20 

10 Xx + xX + Xx = mele) 


10 18 27 





Napkin laundry equations: 


Napkins sent Napkins sent 


tO Fast 


Day Laundry (2) 


10 


Bae Objective function, then, is: 


©. 10(X 


fmable (4.7). 
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The solution of the problem is shown below in 





es 


TABLE 4.7 


SOLUTION OF THE WARDROOM STEWARD PROBLEM 
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$26.70 
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Percan oe. NOecd enae, ,OmM ener cnlra Gay, ewentcy 
soiled napkins are stored, overnight, and not sent to any 
laundry until the following day. The solution would have 
been the same had the twenty napkins been sent to the slow 
laundry on the third day and had then been stored as clean 
napkins on the sixth day when they were returned. In any 
case twenty napkKins have to be stored in one form or another. 
Also, of course, the storage of napkins on the eighth, 
momch and tenth days is in preparation for termination of 
the problem. 

Although this problem spans only a period of ten 
days it could be extended for whatever period of time one 
desires (assuming, of course, one has computational capa- 


city to handle extended problems). 


(Wardroom Steward Problem) 
(Modification No. 1) 


Assume, now, that the Chief Steward finds there is 
another laundry which, for 4¢ each, will launder napkins 
and return them in one day. The Chief Steward wonders if 
he can improve his previous program by taking advantage of 


this new laundry's service. 


| a 
7 





abe: 
To solve this new problem we need only to add an 
additional range of variables (say: w) describing the new 


service; this may be done as follows: 


Napkin usage equations: 


Napkins Napkins Napkins 
Napkins Returned Returned Returned Napkins 

Day Purchased (Fast) (Slow) (Fastest) Required 
(x) + (2) + (y) + (w) = (R) 
1 xe = 50 
2 x5 + Xa = 60 
3 X3 + X43 + X40 = 80 
4 Ky + X19 + X51 + X23 = rae) 
5 Xe + X13 + X55 + Xana = 50 
10 x at x + x a Xx = 100 


L@ 18 Dey 49 





t 
Napkin laundry equations: 


Napkins sent Napkins sent Napkins sent 


to Fast to Slow to Fastest 

Day Laundry Laundry Laundry 
(z) (y) + (w) + 
1 X34 + X54 + Xai + 
2 X15 + X55 + X05 + 
3 X13 +- X53 + Xyn3 + 
4 Xia + X54 + Xana + 
5 Xis + X56 + Xac + 
10 X56 = X36 a Xeo + 


The new objective function is, then: 


0.10(X, +e. -+X1 9) fe 0.03(X,, +e. oe et X54) 


+ 0.04(X +e. . ot X.) 


41 30 


The solution of this new problem is 


table (4.8), below: 
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TABLE 4.8 
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SOLUTION OF THE WARDROOM STEWARD PROBLEM 
(MODIFICATION NO. 
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It can be seen that in this new problem only once 
does the Chief Steward use the one day service, and by 
doing so he saves 20¢ from the previous program. 

It should be noted that, as in extending the prob- 
lem in time, it could just as well be extended also in 
laundry service sources by simply adding more variables. 
The limiting factor would be, as before, the computational 
capacity of the solution equipment being used. 

Since this problem has been cast in a somewhat facet- 
ious setting, one might question the validity of using 
linear programming and a computer to determine how best to 
send napkins to a laundry and thereby save a few pennies. 
It bears repeating that the problem just solved has much 
more serious implications, as mentioned previously (such as 
in repair and maintenance programs). Where we are saving 
only pennies here, by simply displacing decimal points of 
the costs associated with the several variables, it may be 
that millions of dollars could be saved, 

Finally, one other point also bears repeating here. 
AS pointed out in chapter II, the context of a problem is 
not so very important; the structure of the mathematical 
model is, however, all important. It has been shown how 
problem settings may change and the mathematical model 
remains constant: the point being that one must not become 


Steeped and imobilized by classical sterotypes. 





CHAPTER V 
SPECIAL FEATURES OF LINEAR PROGRAMMING 
General 


There are two special features of linear programming 
methods which add to the usefulness and versatility of a 
linear programming approach to problems. These two BARES 
of the technique are the use of linear programming in 
(a) sensitivity analysis, and (b) parametric programming 
studies. 

It has been mentioned that one aspect of data, 
generated within a business organization, is that the data 
are not always completely accurate, that some data must be 
estimated, and that some processes must be assumed. 

A mathematical model cannot discriminate among data 
which are accurate to ten significant figures (for example) 
and data which is not necessarily accurate to even one 
Significant figure. The solution of a model, therefore, 
will naturally reflect only the data used, and in the 
numerical examples given thus far it has been assumed, 
although not specified, that data used was exactly accurate. 

By the use of sensitivity analysis and parametric 
programming one may knowingly use inaccurate data, when 
necessary, and explore the reactions of a mathematical model 


under various assumptions. 
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It must be admitted that the two terms: “sensitiv- 
ity analysis," and "parametric programming," are almost 
synonymous terms and there is some confusion in the litera- 


ture regarding the actual scope intended in this area. 


Sensitivity Analysis 


Garvin is the only authority, which has been found, 
to treat the particular subject of "sensitivity analysis.“~> 
If this concept is considered at all, other authorities 
generally include it in the subject of “parametric pro- 
gramming." 

Garvin?° describes sensitivity analysis as being 
applicable, after an optimum solution of a mathematical 
model has been achieved, and that the concept then encom- 
passes the broad area of five categories, as follows: 

a. Changes in the right hand side of the constrain- 
ing equations (changes in b,, of the general linear pro- 
gramming statement). This includes variations of the 
restrictions on various, or specific, resources. 

b. Changes in the objective function (changes in c,). 


These changes include variations in costs or profit margins. 





>°garvin. Ope 1Cl Geass iGha peer ea. 


-Srbid., p. 49. 
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c. Changes in the variable coefficients (changes in 
a, j)- These changes are relative to the variation in usage 
©®£ resources. 

d. Addition of new variables. This includes 
ehmanges to the product mix and addition or deletion of 
alternative choices in the decision making process. 

e. Addition of constraining equations. This includes 
variations of product dependency on various resources or 
additional usage of different resources. 

Garvin presents the theoretical approach to sensi- 
tivity analysis and how it may be accomplished in hand 
Ponpucations. >. 

The five phase computer program (Linear Program- 
ming II [10.1.008]), used in the solutions already pre- 
sented, may be utilized to accomplish (a) and (b) above, 
(i.e., changes in the right hand side of the constraining 
equations, and changes in the objective function costs, 
or profit). In order to accomplish (c), (da), or (e) above, 
using the computer program, one must solve an entirely new 
problem from the beginning; and in fact, this was done in 
chapter IV (see the "Distribution" and “Wardroom Steward" 


problems). 





Se idae oo. SOne 
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Parametric Programming 


Parametric linear programming may be thought of as 
being part of sensitivity analysis and may be defined as 
the treatment of some element of the mathematical model 
(other than X,) as a parameter of the model. The element, 
now a parameter, is allowed to vary continuously through 
some range of values, while the reaction of the system is 
observed (thus, a parametric study of the system is made). 

Garvin limits his definition of parametric linear 
programming to only a treatment of the right hand side of 
the constraining equations (b, ) as a DeVaeiteere 
Beaty and Gea discuss only changes in the cost 
coefficients (c. ,) of the objective function. It would 
seem that the difference among the authorities, as to what 


parametric programming really is, does not warrant argu- 


ment. 


Example of Sensitivity Analysis and Parametric Programming 


PeOgranmiig , 1 a StracefCerim lone. means wup lanham. 
and, actually, linear programming finds its most useful 


application as a planning device. This idea will be 
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is 
extended somewhat in the following chapter in connection 
with “long-range planning,":; however, now, consider a 


simple example of a planning problem. 
(Planning Model No. 1) 


AnTOLl “company ~witch produces. crude. o1l. in twe 
fields and has a system of pipelines, refineries and 
sales offices to distribute its products to various cus- 
tomers, might have a system as shown in figure (5.1). 

Various company officials have submitted data, as 
their estimates for the following year's operational 
capabilities and prospects; these data are shown in 
Plgure (5.2). 

After constructing a proper mathematical model for 
the data, as originally submitted, and solving the first 
problem, an optimal solution is obtained as shown in 
table (5.1A). 

The planning committee, being in doubt about the 
accuracy of the operational costs of refinery "A", decide 
to investigate the reaction of the mathematical planning 
model by increasing the operational costs of the refinery, 
in various increments. It is discovered that the optimal 
solution remains fixed, although income changes (as shown 
in figure [5.3]), until the operational costs of refinery 


"A" have increased a total of 16¢ per barrel (i.e., froma 
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Aig soa 
SOLUTION OF PLANNING MODEL NO. l 
PL eS 
| | A B 
Description Variable | -Optimum Optimum 
QOuancvit Oliciehesle 
Sauce production from oil | 
mr1eld No, 
iu BE 887.0 3837.0 
2 | 2 0.0 | 
S. 3 EOPOOO HO Ak = OOO. © 
A 4 00040) 18.000. 
5 | 5 7,500.0 | 7,500.0 
Transportation in pipe- 
line No, 
il 6 LOVE S77 JOC RO. 65 76 
2 Bee es isvelo) alee Ls aiielo . 
Refinery usage of pro- 
Qucedad crude 
ferude ) (refinery) | 
(light) (A) 18 4,887.0 | 2,677.0 
(heavy) A) | 10 PR ccylcyeyoyn | Meal AGL 0 
| Gige des B | 9 OL | 2.2 VOe® 
(heav B dele aka fe Se ee) OF Os © 
Crude sales | : 
(ae ht ) oe lee 6000.0 S000. © 
(heavy ) | 13 


Sales of refined products 
(refinery) 





















4,000.0 4,000.0 


OW ete oe ee 


Bld SO oO SO 


4,548.0 OG foes 0 
4,210.0 34-000:20 


kerosene ) (exe 
kerosene (B) : sa pa ey CeO 
fuel oil A ) | 18 4,694.0 2,484.0 

| 19 ls S062 Tap 6. © 
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POcol ret S12 OOD ole acont 2 16/oble) = When ents letter 
GOncitrvon occurs the optimum solution changes as shown in 
table (5.1B). These findings are reasonable since, at 
$1.16/bbl. operating cost in refinery "A", it is less 
expensive to allow refinery "B" (with an operating cost of 
mln la/pol jo to. produce More products. Thererore, a shirt 
of crude oil usage is made between refinery "A" and 
refinery "B" with resulting changes in their respective 
production of refined products. An exact balance is main-— 
tained, however; the crude taken away from refinery "A" 
exactly equals the increased crude given to refinery "B", 
ama tae production lost ac refinery "A" is made up at 
refinery "B". 

It can be noted that the market forecasts indicate 
Peco Emre Orisa s 000 bbls. 8 Gektncd, pucoenees eam be 
sOld: however, it is found to be more profitable not to 
attempt to meet all of the market demands for gasoline 
(even though the revenue is greatest from gasoline.) 
Instead, the optimum solution calls for letting oil field 
No. 2 lie fallow, producing nothing, thereby leaving 4,000 
bbols. of crude raw material unused. The reason for this, 
of course, is that the demands for kerosene and fuel oil 
have been met (with the optimal solution); to produce gaso- 
line in sufficient quantities to meét its market demand 


would cause kerosene and fuel oil, which are produced as 








29 
"by products" in the refining process, to be unsold. It 
ie Wiser icMen. MOt to Use the cruce ol of field Noe 2, 
than it would be not to sell some Kerosene and fuel oil. 
At this point the planning committee might institute further 
parametric studies to determine if additional costs should 
be incurred to develop additional sales for kerosene and 
fuel oil, in order that all available resources could be. 
profitably used. 

The cost changes indicated above were accomplished 
by the use of the “cost changer phase" of the computer pro- 
gram and were made without reworking a completely new 
problem. Similarly, changes to the right hand side of the 
constraining equations may be accomplished by means of the 
"sigmhe hand <sice changer phase” of che computer .predgram, 

By using techniques such as illustrated above a 
business enterprise may discover imoortant reactions of 
its system, which otherwise would not be obvious. 

The planning model, sensitivity analysis and para- 
metric linear programming are not limited to business enter- 
prises, obviously. These methods could, just as well, be 
applied to a supply or logistics system (large or small). 


For example, the following analogy could be drawn: 
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The objective to be optimized, in a logistics system, might 
be the minimizing of costs; however, one could, also, 
optimize (maximize) the number of requisitions successfully 
filled. After reaching one original optimum solution, the 
techniques described in this chapter could be utilized to 
study various reactions that a particular logistics system 


might have under varying system changes. 
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CHAPTER VI 
LINEAR PROGRAMMING AND LONG-RANGE PLANNING 
zt 
General 


The use of linear programming as a planning aid, or 
technique, was touched upon in chapter V. It is the 
intent of this chapter to extend that idea somewhat. 

In many situations, occurring in business, as well 
as in U. S. Navy Supply operations, one is not always 
looking for the one answer to a problem. Generally, what 
is actually being sought, or what is really needed, is a 
range of alternatives, or plans, which can be followed in 
the future when various circumstances are encountered. 

In previous chapters, specific numerical examples of 
linear programming applications have been demonstrated, and 
certainly linear programming has, in recent years, enjoyed 
a rapid rise to fame based on special applications; undoubt- 
edly there will be many other specific uses made of this 
technique in the future. However, it appears that one of 
Ehe Most USerul and Significant arpplleations eo weich 
linear programming may be put, now, and in the future, is 
in an overall treatment of a system as a composite of inter- 


dependent functions. 
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Long-Range Planning 


Although it is not the purpose, here, to fully and 
adequately discuss long-range planning as a subject in 
itself, nevertheless, a short resume of the meaning and 
scope of long-range planning should be profitable, as a 
background. 

Payne describes long-range planning as ". . .one oe 
the really new techniques left to management that gives a 
company a major competitive advantage."°- Wrapp defines 
long-range planning as". . .that activity in a company 
which sets long-term goals for the firm and then proceeds 
to formulate specific plans for attaing these Becika 1Oe 

To produce a long-range plan every facet of a busi- 
ness must be explored (insofar as possible) for five, ten, 
or fifteen years in the future, or for whatever length of 
time is necessary to carry out company objectives. The 
long-range plan in final form will tell management: 

(1) what the company is going to do, (2) how the company 
should proceed to accomplish the objectives, and (3) when 


the company should take various actions to accomplish the 


Objectives. 


Olsayne, Bruce, "Steps in Long-Range Planning," 
Harvard Business Review, (March - April, 1957), p. 95. 


Oe urapp, Edward H., “Organization for Long-Range 
Planning," Harvard Business Review, (January - February, 
LIST iD. oo. 
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There are many aspects to long-range planning which 
must be taken into consideration by an organization desir- 
ing to establish an effective plan. Some of the important 
facets of long-range planning include: 

a anagemene ypebtGy airecki on wm [Op managemene has 
the responsibility and cannot disregard its obligation for 
the establishment of company goals, objectives and policies. 
Not only in this regard but also in all phases of long- 
range planning top management support is essential for 
success. 

b. Personnel orientation. Because long-range plan- 
ning must be comprehensive and tends to “look under all 
corners of the rug" all company personnel must be educated 
to understand the reasons for planning in order to prevent 
pre judicial objections. For a finalized plan to be created 
and become effective, long-range planning must be a “team" 
effort and all personnel must understand (insofar as certain 
sensitivity restrictions or the confidential nature of plans 
will allow): 

(1) The plan itself, and what it is supposed 
tO, accOmoLish. 

(2) The “why" of all the steps and phases of 
the plan. 

C... Aneuivsis se) “To design a workaote "ersecei ve  long— 


range plan a complete and thorough analysis must be made of: 
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(1) Present products. 

(2 eer UruUmce sales pOceite talse 

(3) Market trends and marketing variables. 

(4) Sindustyy growth. 

(5) Company strengths and weakness, manpower, 
finance smprOcUiecrioOn Limca tions.. Cle. 

The planner has to learn the facts by identifying: 

(1) The opportunities and true costs of pro- 
ducts. 

(2) The potential contributions of different 
intra-company operations. 

(3) The economically significant cost centers. 
Once the facts are Known the long-range planner has to 
Know how resources must be allocated to attain the results 
and goals anticipated; thus he must determine: 
© (1) How resources are presently allocated. 

(2) How resources should be allocated in the 
Pewee CO, SuppOre aellviEl eS OL -Greavcese ODPOrtunhity. 

(3) What steps are necessary for the company 
to take to go from present circumstances to what ought to 
be the circumstances for fulfiilment of the anticipated 
goals. 

d. Timing. Once a long-range plan has been designed 


it must be put into a time frame which is: 


a 
_. 
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(1) A realistic framework within which the 
company can effectively work. 
(2) A flexible framework which can be modified 


as exigencies occur. 


This brief outline of the ramifications of long- 
range planning is by no means complete but it should serve 
EO, enOw the extent to which the elements of an organiZza-— 
tion must be examined: in short, a company must search 
out and analyze all variables pertinent to its future 
Sx Seence . 

Planning Or “COur Se 41S Obs Unique COvconmpecLEive 
business enterprises and, as analogies were drawn in the 
previous chapter, similar comparisons can be made between 
commerical operations and military supply endeavors. 

By examining the brief outline of long-range plan- 
ning, given above, it should be evident that the real 
essence of making a long-range plan lies in proper and 
effective analysis of an organization. Without effective 
quantitative analysis a long-range plan (if it could be 
realistically called a plan) would emerge as merely a 
qualitative description of the futuristic hopes of the 


On ganizavion. 
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Lite ina eron Of Planning ane” Programming 


From the perspective of long-range planning and 
business concepts generally, it is a recognized fact that 
the future existence of an organization involves the ful-— 
fillment of specific or multiple goals (which may not always 
bee ee’ Maximizing OL proLits Or MinamiZing Of costs): 
white? presents an interesting essay of possible goals 
which a firm may have, but, as White states, "froma 
Dhaceical “standeoine Goals Of LIirms Or Of “individuals 
typically are expressed in optimal Coenen It is not.so 
widely recognized, however, that an organization must 
"optimize" its operations over the whole of the business 
system. 

The concept of organizational optimization has been 
derived from the simple fact that different alternatives 
to a business decision generally (and usually) involve 
Gifferent costs and/or different variables. It is not enough 
that the quantitative estimates of alternative decisions 
be balanced one against the other, but the impact of each 
alternative must be evaluated within the whole of the 


°Spoulding, Et al. Op. Cit., pp. 181-201. 


C4 tbid., pp. 186. 
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system to produce effective long-range plans. The danger 
to be avoided is "“suboptimization," which results when 
alternatives are considered independently of the system. 
Suboptimization is a case of optimization for one phase 
(or element) of an operation, without taking into considera- 
tion every factor which has a bearing on the problem, 
whether in an obvious or in a subtle way. 

TO 7Overcome “ene pDierall Gf —subopecimizaticncvelis 
necessary to look simultaneously at all pertinent elements 
of a system as they progress through time and affect one 
another. Within its limitations, linear programming offers 
an ideal approach to long-range planning, while overcoming 
SUbODeImaZation. 

To illustrate what may (or must) be done in the 
realm of long-range planning, the example (planning model 
No. 1) introduced in chapter V, will be extended. (Fig- 
ures 6.1 and 6.2 and tasle 6.1A may be thought of as 
extensions of figures 5.1 and 5.2 and table 5.1A respect- 


ively.) 
(Planning Model No. 2) 


Assume that the planning comnittee of the hypo- 
thetical company, of the earlier example, wish to look 
ahead for not only one year but for two years. Due to 


various contingencies, the company will have accomplished 
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certain modifications in its operations. The second year's 
Situation is shown in figure 6.1, and the prospective opera-— 
tional and investment data are shown in figure 6.2. 

It 1s convenient when discussing linear programming 
planning models to distinguish between two types of vari- 
ables: (1) operational variables and (2) investment vari- 
ables. 

Operational variables are associated with the 
various operational elements of a business system, and 
usually are concerned only with the cost of operating the 
particular business. Investment variables are generally 
related to the industrial plant and facilities and pertain 
to capital expenditures. Operational variables describe 
the operations of a business for a relatively short period 
of time. Investment variables link operations through time 
periods (or within a time period) by describing the effect 
of new and/or different equipment on operations. 

In the previous example (planning model No. 1) only 
operational variables (x) were used. In the present 
example one may notice, in addition to the operational 
Var toaples)."he introduction cf investment. variables (Y 5). 
By using investment variables the planning model is given 
what might be called, and in fact is a new "dimension". 

The investment variables are not limited to represent-— 


ing only realities; actually their real usefulness lies in 
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TABLE 6.1 


SOLUTION OF PLANNING MODEL NO. 2 

















| 





| 7 a = 2 
Description po as Optimum pptiman 
lcuantity Quanti Veni ey 










Crude Production [rom 
oil field No 


ili ees Spas) OS en ©) 
aa 
| wy i Ree tel 4,500. 4,500.0 | 
3, O00. O | 5 ae 81 O16 eG 


O7OD es om 7, 200 20 






833.0] 3,500.0 
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oreo) O10 58 









Prensportcecion. an 
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38 | | 
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Increment to pipeline #2 |} 
Refinery usage of pro- 
duced crude in 
refinery "A" 
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Crude sales 
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Increment to sales of 
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| Sales of refined PLOcuEGES 
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representing proposed system modifications which may, or may 
not come into being in the future. 

In the present example, the planning committee has 


introduced two investment variables (Y Yoo) to determine 


49! 
the advisability of increasing the cdlemands for certain 
products. The investment variables represent the elements 
necessary for increasing demands, such as: hiring and 
training of additional salesmen, providing more advertis-~ 
Hage. Miiereasing Product. appeals erecwyiin teble Go, 1A Atrcan 
be seen that the optimal solution calls for expenditures to 
increase demands for refined fuel oil (Yeo) Jeon eam gleqsuum ese 
increasing the demands for light crude oil (Yy9)- 


One investment variable (Y has been weed to 


43? 
investigate increasing the capacity of the pipeline facili- 


ties, and one investment variable (Y,.) has been used to 


47 
indicate a possible increase in the processing capacity 
Bester rcevOn-Or rebrinery A 

The simplest use for the pianning model is deter- 
Ming hee SO liGIOneInvVel ing. Pies GONG elem s (suchas. lee me 
previous example (planning model No. 1) where the solu- 
BLOM: May noe raCeeCEea1dce eo ING exoCe wre tie sCOoperon ene 
planning model, to be really efrfrective, must be extended 
beyond these limits since conditions are seldom if ever 
SEAELC> NOWEVer ea WOE OL “Caution 1S an,order. “By the 
introduction of fictitious elements, or elements which 


later may be discarded, one has automatically given the 
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POMS on cibwea Geer ta iis ScoMIIO US Mea cumcn sla ©lNer PwOrcs,.. in 
these situations one cannot, necessarily, accept the 
"answer" as being literal: the "answer" must be examined 
Mien regard co the conditions imposed on the model. In 
Pies FeSpec. Sensitivity analysis and parametric program 
ming may be used to great advantage. For example, the 


planning committee, using investment variable (Y,—) made. 


47 
Mice citer al. assumption chat, above 14,000 bbls che 
Capacity orc refinery "A" could be increased at a cost of 
BOC eC er ODl, *wby Increasing CcOSEs, €O 25¢ > Der polleic ne 
found (as shown in table 6.1B) that the optimal solution 
changes not only for the investment variable (Y 7) but 
also, among other changes, for investment variable (Yoo) 
representing proposed expenditures for increasing demands 
for fuel oil. Similerly, by allowing demands for gasoline 
te cecrease from 16,000 bbls. to 8,000 bbls.,, it is found 
that there is no change in the optimal solution until 
demands for gasoline decline to 8,000 bbls. (as shown in 
table 6.1C). At this latter point the new optimal solu- 
tion calls for less capacity of tne refinery and less 
expenditure for the development of new demands for refined 
Fuel oil. The interaction and interdependency of the 
variables, of course, account for changes among the vari- 
ables wnich are only indirectly related. The changes in 


variables (Y 4) and (X as shown above, prevent 


aa). 
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“suboptimization" because the system, rather than individual 
elements of the system, have been optimized. In addition, 
it becomes obvious that any one solution is devendent upon 
many factors and the static conditions prevelent at the 
time; given some change in fixed conditions the solution 
may, and probably will, be different. 

in addition to investigation of real and hypotheti—- 
cal elements of a system, as indicated above, the planning 
model may be used for many other planning analyses. How- 
ever, the planning model cannot describe anything about a 
great many things such as the criteria on which assump- 
tions were based in ord2r to formulate the model; these 
criteria must be set by management using experience and 
judgement. Also, there are a great number of other factors 
involved in business which are (at present, anyway) non- 
QMiaveiriabte. Such Lactors aneclvces cenpany “qooowl 1. 
labor relations, public relations, ethical practices, 
legalities, etc. As Rapoport and Drews have stated "the 
establishment of criteria as to what is actually best for 
eeouciness (enter orice we a btogetnes beyond —s= weaimacr 
mathematics,"°> 

Although in the example usea, here, the time span 
covered has been only two years, it could have been extended 


©°R apoport and. Drews. Oo. ©2the. oO. 79. 
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Well pebevyora het. “fae Jimicing fLaccor 1s. ageim as we have 
seen before, the computational capacity of the equipment 
eing used, For example, Rapoport and Drews estimare that 
a realistic planning model, covering five successive time 
pert tods) would contain, ~verhaps 500 verniables and aoout 


66 : he : , 
i A moagest model oF that size would require 


400 equations. 
considerable computational cépacity, but even with this and 
other disadvantages, a .inear programming approach to long— 
range planning offers decision makers a very powerful tool 
POr onmalysis in. determining PulLure courses action. 


Oe be wo. 


Cie eis “Vw 


SoM UNVESTICATIONS iro EXTENSIONS 


OF LINEAR PROGRAMMING 
General 


During research for this thesis an attempt has been 
mace Fy Liaienauunor . Ode VoLoD co a limleed extent, Ga 
tain new and different approaches to two linear program- 
ming problems. This chapter consists of the results of 


the investigations in this area. 
me WOagrng Probie 


The loading problem, whicn is of some military value 
if not also of some theoretical interest, has been 
MetietOneA yoy OnemMeuchor es Neving Nad nO ws 1 oni leant jeony 
geveancea LOr Les) SsOlLUrIOn jas CF hereon 

The loading prob.em may be szated simply as follows: 
Suppose some vehicle (an airplane oz ship, for example) is 
tO be Loaded up to or including &@ certain maximum weight 
(W). The cargo which is to be load2d is in the form of 
discrete packages, or bundies, eacn having individual 
values (v.) and individual weights (w,). The problem pre- 
sented, then, is how to load the vehicle with the maximum 


: o7 ;, . 
: Gale, David. Tne Theory of Linear Economic Models. 


New York: McGraw-Hill 300k Co., Inc., 1960. p. 134, 
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value of cargo but not exceeding the maximum allowable 
weight load limit (W). 

The problem has been stated in linear programming 
68 


ff OCM: 


Maximize: 


ay Ez Ve 
ies I ~ 
oe eee oe ee J CO 
> L 
SUD JeCl reo: n 
7 be = 
ee iMG ys 


Present known methods of solving linear program- 
ming problems, such as the simplex method, cannot be used 
to solve this problem since they do not guarantee integer 
results. Even the transportation method of solution, which 
yields integer values to certain tyoes of linear program- 
ming problems where the problem may be formulated by means 
of double description or double identification of tre 
variables, cannot be used. The transportation method of 
solution requires known source availability values and 
destination requirement values as well as the various 
costs associated with moving an item from one location to 
another. In the loading problem variables cannot be 
Gescribed in analogous terms suitable for submission to 
eee Eransporcaciom alqoritam. 


OS TDi ds 





148 

The loading problem may be solved, in one way, by 
direct enumeration. Here one merely examines, in turn, 
each of the possible sets of objects to be loaded, and 
among those which do not exceed the weight limit, one 
chooses those which will yield the greatest total value 
for the combined cargo. 

The CGuperculey involved. in aceempting £o compare 
all possible combinations is the magnitude of the number 
of comparisons required. The total number of combinations 
Gin LTMIings taken: (Zero ati a time, one ata time, Ewo 
ae folie eels Ce OUC hn atria tinea us De Cor Dye Is 
if zero at a time is not considered). It is easily seen 
that the total number of possible combinations rapidly 
becomes very large and quickly outgrows our time and 
ability for making comparisons (for example, just ten 
hackages would require 1,024 comparisons to assure finding 
the optimum loading schedule; 30 packages would require 
1,073,741,824 comparisons, etc.). The direct enumeration 
method, therefore, is not a practical method of solution 
for the loading problem (although direct enumeration will 
be used to check the procedure advanced below). 

The author now presents the Following original 
heuristic development as a procedure for solution of the 


loading problem. 
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Be re-examination of the problem statement it can be 


seen that the restriction 


May be restated, after all E. are set equal to unity, as 
a 


follows: 


hus 
= 
! 
” 
I 
= 


i 


Where "S" is some slack amount which cannot be loaded. The 
problem, then, could be restated by saying that we desire 

to minimize "“S". We have, as yet, said nothing about the 
values (v.) of the individual weights (w. ) and we must 

have the two variables related in some way. If we introduce 
an additional variable "X", which shall be defined: 


A 
1 


v/w. Pee ee le) 


Ww. 
a. 


Vi/X, Boe. ee eee ee 
and which will be called, simply, an index, we will then 
have a connection between weights and values. 


Using equation (7.2) the problem may, again, be restated 


as follows: 


Minimize: S 
n 
where: oy v./X. - S =W 
; nts ale 
le 
n 
oa v./X, 
0 
n 
S = W+ XY v [X; 





IFS 18, 
Then, for "S" to be a minimum the derivative of "S" with 
respect to each individual value (v.), (ds/dv, ), must 


approach zero. Taking the derivative we obtain: 


n 
a(-W + «2 v,/X, ) 


as/av . = 
dv. 
i 
n 
NN By oe) 
then: as/avenae=) a (SES ar 
~ ay! 
i) 
eee 
or: as/dv . = av, (v,/X, tV5/Xot+V3/X5+- 3 tv ./X, ) 
then: as/dv ; = 1/x, 


For the derivative (dS/dv, ) to approach zero, each individ- 
ual index (X, ) must be as large as possible. In other 
words, by choosing bundles having the largest indexes we 
should, theoretically, have a loading schedule containing 
the greatest value of cargo with the minimum amount of 
weight. 

If, then, we are given a listing of various weights 
and their associated values, by dividing each value by its 
respective weight we may obtain an index. Arranging the 
various indexes in descending order of magnitude ina list- 
ing with the respective weights and values we may choose, 
from the beginning of the list, those weights which add to, 
or are just less than the desired maximum load limit. For 


example, if we were given the following weights and values: 





a new listing, arranged by descending index, 


formed, as follows: 


PACKAGE 
(a) 
(b) 
Kon) 


Then, if we were given various loading objectives, 


We tGHT. VALUE 


OR LOe 
iM i 
a iG? 
Ze ee 


WE LGHT VALUE 


Za IE Oes 
De Ge 
ine is 
LOE LOX 


alpoys 


loading 


schedules containing the optimum load could be prepared, as 


follows: 


MAXIMUM WEIGHT 
(LOAD OBJECTIVE) 


2 
S 
4 
2) 
fi 


10 


LOADING SCHEDULE 


(PACKAGES ) 
(a) 
(20) Fey, 
(aie Ce) 
(20 ae Ce, 
Catan 
G2U) me Oo) yan ey), 


LOADING WEIGHT 


(OPTIMUM WT, ) 


2 


CO iN [Ww [WwW Iw 


i. 
20. 
Ze 











It should be noted that where indexes are identical 


the smalles weight should precede the larger weight in the 


list (as between packages (c) and (d)); where a weight 


exceeds the remaining allowable load it must be skipped and 





da 
the remaining weights chosen which are equal to or less than 
the remaining allowable load (as shown in obtaining the 
loading schedule for weight objectives 3, 4, and 5). 

At this point one has obtained a loading schedule 
which is equal to or less than the maximum allowable load 
limit and which has what might be called the highest "value 
density" of all possible loads. However, since we must 
find the highest absolute value of the cargo within the 
allowable range, the load having the highest "value density" 
will not necessarily be the loading schedule having the 
highest total absolute value. This may be explained by 
referring to the initial discussion concerning the mini- 
mizing of the slack (S). By taking the derivative 
(as/dv, ) we have tacitly understood that we have an infin- 
ite number of bundles, or packages, of infinitesimally 
small incremental weight and values. Such, unfortunately, 
is not the case, as a practical matter; if it were we would 
have already solved the loading problem. 

Once, having obtained the load of highest "value 
density" it may be possible to improve on this loading 
schedule by replacing some weight or weights on the schedule 
by some weight or weights not included in the loading 
schedule. 

The algorithm may be followed, most easily, ina 


step-by-step process, as follows: 


- 
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step 1. An index (X,) is computed for all data 
POInes. whieh Consist, or various weights (w. ) and associ- 
ated values OUR BIE 


X, = vy/w, 


step 2. All data are then sorted into a list by 
descending order of index. 

Sten 6. - ihe weighe (w, ) of the first data point is 
compared with the maximum allowable load limit (WMAX). 
There are three possible courses of action: 

a. If: w. < WMAX, then the data point (wi, ee is 
put into a list (to be called the "loading schedule"). 

One would then go to step 4. 

Pere “Lie Ww, > WMAX, then the data point (w. , v,) 
is put into a list (to be called the "excluded list"). 
One then would go to step 4, 

oe oS w, = WMAX, then the data point (w., vy) 
is placed in the loading schedule and all other remaining 
data points are placed in the excluded list in descending 
order of index. One then would go to step 5. 

Step 4. The quantity (WMAX) is reduced by the amount 
of the weight (w, ) just placed in the loading schedule, so 
that WMAX now equals the balance which can be loaded. One 
now repeats step 3, using, in order, the second, third, 


fourth, etc., data points until all data have been exhausted 


in the comparison with WMAX. 
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At this point one has obtained a first feasible solu- 
tion. A checking procedure is next carried out to deter- 
mine if this first feasible solution may be improved. 

Step 5. The data points (we, v,) of the loading 
schedule are sorted into a listing which is in order by 
ascending values, (v.), where v, is the smallest value on 
the loading schedule and v, is the largest value). 

For simplicity we will denote the data points on 
the loading schedule by the additional subscript (s) such 


as (we Vv 4? and we will denote the data points on 
? 


_ @ 
i Ss, 


the excluded list by the additional subscript (0), such 


as (w ). If, then, "n" equals the number of data 


eZ ; 
O, J O, J 


points on the loading schedule, and "N" equals the number 
Of cata points on the excluded lict )thenes (i-1,2,4 02. n) 
EC j= Zoe 6 eG) 


Step 6. The value Oe ;? where j=l for the first 


f 


data point on the excluded list is compared with the value 
eae i? where i=l for the first data point on the loading 


senedule, vand Jtetting: Vv . = VS, and w . + S = WS, 
S, l Ss, l 


where "S" eguals the difference between the maximum allow- 
able load limit and the total weight on the loading 


schedule, then the alternatives are: 


eee aa = VS, then no replacement can be 


oO, J 


made, and one would go to step 9. 





USS 


Deen ies ve F > VS, a replacement may be possible, 


and one would go to step 7. 


Step 7. 
et 
ae Lf: we ue wS, then vey a) and Wey a! are 


placed at the end of the loading schedule (the n+l position). 
After the data point is placed at the end of the loading 
schedule, a check is made to determine if any of the data 
points which make up WS can remain on the loading schedule, 
Then, those data points included in WS which have to be 
removed from the loading schedule are set equal to zero. 
One then would go to step 9. 

Dba 2br saw . > WS, then we make further compari- 


O, J 


sons; one would go to step 8. 


Scope G. 


a. If the value of the next data point wz 


) 


bade 


is less than the value of the excluded point being con- 


sidered: 


— ., then one would add the data 
s, itl OR ks 


Vv 


eee s41)? on the loading schedule, to WS 
a 


eleueue (Ww. 


é 
and VS respectively, so thats 


WS ew = WSold ‘ “s, i+1 


VS new rs Vola a ‘s, i+l 


Then, one would follow either step 6a or step 6b, as 


appropriate. 


| 
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IV 


V then one would take the 


lope! Geen 4 ; = 
s, itl ~ Yo, 3 


value of the next data point and repeat step 8. 

Go 1&2) °Step 6 1S reached and all data points on 
the loading schedule have been considered in comparison 
with the data point OG, j, ve y? on the excluded list, 


such that we j > WS, then no exchange between the two 
é 
lists is made and one would go to step 9. 


V ya 


Step 9. The next data point SOT re See 


the excluded list is then used in the same way as point 
Vora von and one then would follow step 6 from the 
beginning. This procedure is followed until all of the 
data in the excluded list have been checked. 

step 10. After the above procedure is followed, 
the answer has been obtained, as shown by the loading 
schedule. 

A computer program (written in Fortran II language) 
was developed in an attempt to solve the loading problem, 
as presented above (for simplicity this program will be 
referred to as the approximation algorithm). In addition 
four other computer programs were prepared in con junction 
with attempting to determine the accuracy of the approxi- 
mation algorithm. These latter four programs consist of: 

a. A random data generator which prepares simulated 


weight and value data and simulated maximum load limits. 





ies 7) 

b. A direct enumeration program which makes all 
possible comparisons for six weights and six values (64 
computations). 

c. A direct enumeration program which makes all 
possible comparisons for ten weights and ten values 
(1,024 computations). 

dad. A direct enumeration program which makes all 
possible comparisons for twelve weights and twelve values 
(4,096 computations). 

Copies of all the above mentioned programs may be found 
in Appendix I. 

The simulated data, prepared from the random data 
generator program, was submitted first to the appropriate 
direct enumeration program to determine the correct 
answer and then to the approximation algorithm to deter- 
mine its accuracy. The results of these tests are shown 


below in table 7.1 and table 7.2. 


TABLE 7 oi 


REOULTS OF TRoTS CONDUCTAD WITH APPROXIMATION ALGORTIN 


| INCORRECT 
NUMBER OF | TOTAL NUMBER | ANSWERS FROM 
| DATA POINTS | OF TEST SETS | APPROXIMATION | PER CENT ERROR 
| ALGORITHM 
6 38 4 | OsS 


12 1 10.0 





eos 
Dice EE eee 


ANALYSIS OF ERRORS MADE BY APPROXIMATION 
DOGORT IRM IN EIGHTY —-ON® ESTs 


Weight Results 


Weight by 


ae Maximum apne a Approxima-j Weight 
emacs Pe Roose se Difference | Per Cent 
Points Weight 
-128.0 -12.8 
679.0 573.0 lowe 2710-5 | 
= LO 919.1 8383.0 67 10 -167.0 -20.0 
DOg.2 BG eee 904.0 1 - 29.0 - 5.4 
ACM CO CCR -313.0 -18.8 
6| 10 - 5.8 
eee’ 407.6} 392.0 |! 292.0 -100.0 -25.5 
a eee i ee ee es merce 
9 2s _}i.as.s} 2075.0 12 Vow eb eeee oF? ee 
6 _[a,437,6[1,227.0 | 904. 
fae gt 750.0 | 342.0 | -408.0 54.5 
Value Resuits 
| | Number St am neh pl ae EE | 
No. eee tion tion, Algo-| Pit tere ices! Perm ecans 
(Points| Value ngatim)gyic 
r[ 22 | | 2,932.0 
| 10 1,646.0 |1,457.0 | -189 - 8.7 
2 1Q E98 -O lL 240.0 - 41 
2,396.0 |2,378.0 | - 18 2 69s 
5 10 POL B17 Ol Ae (Oe aes 
ne ee ee eee 
el 6 |  }2,200.0 |1,645.6 |-555 | -25.2 
9 6 2- O35 -6 12702020 -312 -10.6 
ie 6 DP 2OOLO 2 Zee 0 - 71 - 3.2 
al A 1,118.0 980.0 -138 Si | 





Approximation algorithm weight exceeded optimum 
weignt determined by direct enumeration. 
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Example solutions of five sample problems, together 
with direct enumeration computations, may be found in 
PO mer Cake i. 

It should be pointed out that data which were used in 
the above tests were of a realistic nature: no data were 
included in the tests which would have resulted in "no 
eonrest'. 

Although these tests are not entirely conclusive it 
appears that the approximation algorithm tends toward being 
exactly accurate Approximately 80, Of. the time and approxi4 
mately 20% of the timé it would be within approximately 75% 
of the maximum possible loading value. epending, of course, 


on the accuracy which is required in determining the maxi- 


thes LOading approximation aloqorithm would be suriveientiy 
accurate. The worth of the approximation algorithm is 
enhanced when computing times (utilizing an IBM 1620 com- 


puter) are considered: (see table 7.3, below). 
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TA BED? so 


COMPARATIVE SOLUTION TIMES FOR LOADING PROBLEMS 


——_ ne eat gg tl 





rade. Computation Total time | 
of Total No. | Time per Total Time |by approx. 
D Gr Set by Dir- {by Direct algorithm 
ata 
Data Sets | ect Enumer- jEnumeration jifor total 
Points 
ation (t.) data sets 
Fig ese OF india 2 
10 33 3.5 min. 2 hours 5.5 min, 


12 10 2 hours 2 min. 


Although computing times increase with the approximation 
algorithm, as the number of items being considered 
increases; however, the increase per data point is on the 
©rder Of microsseconds, Om the oeher hand, by dipect 
enumeration, which, of course, will give an absolutely 
correct answer, computing times are slightly more than 
doubled with each additional item, or data point, considered 
(see figure 7.0). 

For example, to compare twelve items the computing 
Eime WS approximately -IS"minutes- io -cOmpare Emir teen teens 
the computing time would be approximately 30 minutes. It 
is interesting to note that, by direct enumeration, with 
constant computation, the time required for comparing only 
thirty-one items would be in excess of fifteen years, using 


ag eM O20 cOMmeucer., an figure 7 50 eee be Seen liar 
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the computation time increase at a somewhat greater rate 
than the number of comparisons being made. 

The loading problem offers many interesting aspects, 
and during research for this thesis, it has been approached 
from several points of view, most of which proved fruitless. 

It is believed, although this author cannot prove 
his belief, that the loading problem has no direct solution 
(other than the direct enumeration method). The author has 
been led to this belief after considering that the data 
with which one deals in the loading problem, must be con- 
sidered as being completely random, and, therefore, any one 
direct method which would solve the problem in one case 
would, most probably, fail in another case when tlre data 
assumed a different distribution pattern. In addition, 
there are an infinite number of combinations (when decimals 
are allowed) which can combine to make some given total, 
and any number (say, X, ) which is in the set of data and 
which is less than the maximum allowable load (say, M) 
belongs to at least two sets of numbers. The number (X, ) 
belongs tc both an infinite and a finite set of numbers 
which exactly totals to the objective (M); but some of the 
other numbers in a common set with (X, ) are not in the 
available data. Tne number (X, ) also belongs to a finite 
and an infinite set of numbers which totals to some value 
less than the objective (M), and the other numbers of this 


common set may or may not be in the available data. 
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The other numbers in the list of data (say, 
KX, X3,- é -X) which is available belong to similar 
sets of numbers as does X)- Then, essentially the loading 
problem is one of finding among all the available data, a 
common finite set of numbers which either sums exactly to 
(M): 


xX, = M 
a 


n 

i1=1 
Or a common finite set of numbers, from the available data, 
which sums exactly to some number less than (M)s: 


xX. =M-S 


a 3 
ome = 


Wt 3 


i 

Therefore, what one has to determine is one finite set of 
numbers from an infinite number cof sets, without knowing 
what set of numbers one is looking for. In essence this is 
analogous to, but worse than, "looking for a needle in a 
naystack''. 

The author believes, and here again he cannot prove 
this, that the loading problem may eventually be solved 
by some trial and error procedure. One procedure, which 
was not attempted, but which would seem to have good 
possibilities for almost complete accuracy, is a method 
utilizing sampling theory. A computer program could be 
developed to compute (within the desired accuracy) the 


number of samples one would require and then proceed to take 
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random samples of the data for the required number of times, 
or for some very large number of times, with the best 
answer, of course, being stored while all others would be 
discarded. It would seem, also, that a sampling approach 


could be applied to other intractable problems as well. 
The Loading Problem Including Cube (Volume) Considerations 


During the course of research and study of the load- 
ing problem the idea occurred that it might be extended to 
include cube (or volume) considerations in addition to 
just merely weight restriction. 

In this new problem, as in the preceding one, imagine 
that a vehicle (such as a ship or airplane) is to be loaded 
up to or including a certain maximum weight (W) but also 
not exceeding a maximum allowable volume, or cube (C). 
Again, the cargo which is to be loaded is in the form of 
discrete packages, or bundles, each having individual 
values (v.), individual weights (w. ) and individual volumes, 
or cube iC, ). 

The problem presented is how to load the vehicle with 
the maximum value of cargo but not exceeding the maximum 
allowable weight load limit (W) and also not exceeding the 
maximum allowable volume or cube limit (C). 

Although not found elsewhere specified as a linear 


programming problem, this problem will, here, be treated in 





Gs 
the same manner as the preceding loading problem which had 
only weight limitations. 

In linear programming form we might present the 


problem as follows: 


n 
Maximize : ay Baas 

ent 

“ea 
sub ject tos EB; = 1 Or E; = 

n 

ny kB , w. = W 

iz. 7? ? 

nN 

ie Bk {C2 Vane 

ion ee 


For the same reasons given before this problem can- 
not be solved by present known methods of solving linear 
programming problems. Of course, as before, and with the 
same obvious drawback of too extensive comparisons, the 
problem might be solved by direct enumeration. 

In attempting a solution of this new problem we 


might begin by restating the restrictions: 
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after all E. are set equal to unity, ass: 


n 
a w. - S = W 
i=l ; 
n 
y Cap So t= 
i=l z 
where "S" and "S'" are slack amounts of weight and cube, 


respectively, which cannot be loaded on the vehicle. Then, 
if we define a new variable "X", which shall be called an 


index, ass: 


X, = v,/(c, )(w, ) ewe does 
w, = vi/(c,)(X,) . 2. - 2. - ~ (7.4) 
oa v,/(w, )(X, ) a a ae ee, 


we have a relationship among the respective weights, values, 
and volumes of the several packages. 

If we approach this problem as one of minimizing 
the slack we may then state the problem, using equations 


(7.4) and (7.5) as follows: 


Minimizes: S$ and Ss! 
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where: 
n 
x vi /(c,)(X, ) -S = WwW 
Aa 
n 
-S = W - a v7 (ec, )(X, ) 
n 
S= We E vi/(o5)(X) 
n 
2 v,/(w, )(X, ) =—S = 7¢ 
eae 
n 
-S' = aa for v,/(w, )(X, ) 
Then for "S" and ugti to both be a minimum the deriva- 


tives of "“S“ and “S'" with respect to the individual values 


(vi), (as/dv, ) and (dS '/dv, ), must each approach zero. 


Taking these derivatives we obtain: 


ds/dv . 


ds/dv., 


ds/av 


dv 


v./(c, )(X, ) 


! 
=, 
-} 

WS 
rae 


(ee 2: 


vi/(c, )(X,) 


vile, )(X,) + vo/(c,)(X,) eee 


+e a v3/ (cz) (X3) aoe 


tee Vy/ (C4) (Ky) 





ies 








ds/dv , = 1/(c, )(X, ) 

, d : 
ds fav, = Aue -C + x vi /(w, )(X, ) 

1 deal 
d 
ds ' /dav, = dv, an Vv /(w, )(X, ) 
aS'/avs = By | Vy /(\)OG) + v2/(w2)(X) + 
+. ot Waa Viele) Xe) 

ds ' /dv = 1/(w, )(X,) 


For the derivatives (dsS/dv, ) and (ds '/dv, ) to approach 
zero individual indexes (X, ) must be as large as possible. 
In other words, in this problem as in the preceding one, 
by choosing bundles having tre largest indexes we should, 
theoretically, have a loading schedule containing the 
greatest value of cargo with the minimum amount of weight 
and the minimum amount of volume. 

Without providing a demonstration of a numerical 
example at this point (numerical examples are given in 
Appendix II) let it be sufficient to point out that in this 
problem, a first feasible solution is obtained in the same 
manner as in the preceding problem. Packages having weights 


and volumes less than any remaining allowable load weight 





IMS se 
and allowable load volume are chosen by index number in 
Gescending sequence. The first feasible solution is the 
load of highest "value density". 

For the same reason given in the preceding loading 
problem the first feasible solution having the highest 
“value density" is not necesarily the solution having the 
maximum possible absolute value. It may be possible to 
improve on this loading schedule by replacing some weight 
Or weights on the schedule by some weight or weights not 
included in the loading schedule. 

The step by step procedure of the loading algorithm 
with volume considerations is the same as for the previous 


algorithm except that the index (X, ) is computed, as follows: 


X, = v./ (w. )(c, ) 


and additional methods of comparison for (c, ) are used 
which are the same as were shown for (w,) in the previous 
algorithm. 

A computer program (written in Fortran II language) 
was developed to solve this loading problem containing 
volume restrictions as presented above. The computer pro- 
gram is included in Appendix I. Five numerical examples 
solved by this algorithm, and by direct enumeration, may 


be found in Appendix II. Although this algorithm has not 
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been tested in the manner in which the loading problem was 
tested, it is believed that similar results would occur as 


with the previous tests. 


Transformation as a Means Toward Optimization 


Unfortunately, not all (or even a majority) of 
relationships encountered in real life situations occur in 
neat linear form. A portion of the research that has been 
conducted in linear programming has been devoted to treat- 
ing some non-linear relationships, by various approximating 
techniques, as linear relationships (for example, see 
Charnes and Cooper®”). A great number of non-linear 
relationships are of such a nature, however, that they may 
be readily transformed, in some manner, into linear rela- 
tionships not requiring approximation. It would seem, 
therefore, that relationships of this nature could be 
treated, upon transformation, by standard linear program- 
ming techniques when optimation is desired, 

Two types of relationships, which are not only very 
useful but also frequently met in engineering and industrial 
enterprises, are those described by exponential and hyper- 


bolic curves. These relationships take the general form: 


EXPONENTIAL CURVE HYPERBOLIC CURVE 
a aoe vox ax? 


69 
Charnes & (COOper ) (Op. Cit, cChapten ox. 
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These curves upon transformation take the following general 


forms: 
EXPONENTIAL CURVE HYPER BOLI CACURVE 
ine ye = ina -— bx In y = ln a-b ln x 
In y + bx = ln a In y + b iIn-x = Inve 


The exponential curve becomes linear in semi-logarithmic 
coordinates and the hyperbolic curve is linear in log- 
arithmic coordinates. 

' In the following two sections simulated maximiza- 
tion problems of the two general types shown above are 
presented utilizing the simplex algorithm to reach optimum 
solutions. Simulated problems involving only two vari- 
ables have been chosen in order to present graphical as 
well as computer solutions. 

Although multi-variable problems have not been con- 
sidered it would appear that, as long as linear trans- 
formation is possible, there should be no prohibition to 
the use of linear programming techniques to optimization of 
multi-variable non-linear problems. For example, there 
does not seem to be any reason why relationships of the 
following general form could not also be optimized by the 


simplex linear programming algorithms: 
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transofrmed to: 
Iny=Ilina-bi1nx+dtIncec-ti1nw+i1ng-v inz 
Uy eer iescendr ete ee Wot oe eee ng Wied ee Glemeet ae inIn a Gy 
where: 
a, b, c, g, t, and v are some known constants 
WwW, X, Y, and z are some unknown variables. 


By lettings 


in y = X5 
UB ete, «xe — Xo 
in’ w = X3 
iy. 22, xX, 


Pina slay C + lia (Gack 
then s 


X) + bX. + tX3 + VX, = kK 


Therefore, in the above manner, the standard linear program- 


ming problem formulation can be achieved. 


Exponential (Semi-Logarithmic) Maximization Problem 


Assume the following abstract systems 


Oe 2 ox 
e 


Objective (Maximize): z= y subject to 


the following restrictions: 





HWA 


8 


A IA 


HA 


3) 


HA 


IA 


6 


IV 


IV 


Ie, 


S35 


220 


-0.4024x 
[= 


-0.384x 
e : 


@ 70 -230x 


-0.3584x 
e es 


-0O,173x 
e e 


gee: 
(7.6.1a) 
(7 Jeraniay) 
(7.6.1c) 
(7.6.1d) 
(7.6.1le) 


(7e6el ie) 


The above system may be transformed by the use of log- 


arithms, 


into a linear system as follows: 


Objective (Maximize): 2Z 


in 
ln 
in 
Ky 
in 
ln 


By allowi 


a 


KK 


ng 


ng: 


~ 


0.462x 
0.384x 
O.1734& 


0.230x 


Oe2358s 


NA TA TA HA 


NA 


Tri 


in 


im 


8 e o 
OS e 


3S oe 


o>) 
© 
0 


< 
II 
a 


Iba 


Zz 


in y + 0.231x 


(7.6.2a) 
(7.6.2b) 
(726226) 
(7.6.2d) 
(756522) 


(762 ey 
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the linear system, as plotted on semi-logarithmic coordin- 
ate Grapnm paper ian Ligure 7.17 is as Follows: 


Objective (Maximize): Z2 = X. + 0,231X 


it 2 
X, + 0.462x, 07944 eee ee Gee) 
X, + 0.384Xx, 2302501 eee eee. 7G) 
X, + 0.173x, Se 2096 teen 7 eGo) 


HA 


wet we ecu sOOGA SPs « “ho S. Sa a ape Meu ba ae Ov oy) 


Z 


IA 


0.69315 t e e e @ e e e e e (7.6. 3e) 


< 


X, + 0.358x E7O1I6 s oso «se ee mtr CSE) 


I Z 
By the introduction of slack variables the inequations are 
converted to equalities suitable for submission to the 


simplex algorithms 


xX + 0.231x, = Z (Maximize) Objective 
Xy + 0.462x. + X3 = 2.07944 
Xy + 0. 384x, + Xy = 2.30259 
Xy + 0.173X, + Xe =e 20890 
x, + 0. 230K. + Xe = 1.60944 
xX) + Xo = Ow Ose 
xy + 0. 358X, + Xo =e oO 


The above data, in the proper computer program format, is 
shown in Appendix II along with the computer solution to 


the proglem., The computer output indicates the solution 
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to this problem to be: 


X, = 0.608945424; X 3.00799060; FUNCTIONAL = 1,38444540, 


1 2 
ReCGaliln hop atmeit Z = eee 
Ay = in y 
x, = x 
a — ae 9 6 eA 
we can see, thens 
In y = 0.68945424 
or y = 268945424 = 1.997 
and x = 3.07 
m > 
Z; = (1.997) (el 0-231 JL 3-00 ns 
Z = (1.99")(1.997) = 4.07 


But since the functional is given as: 1.384444540 we take 
the anti-logrithm and find it to be 3.99. 

Comparing these results with the results achieved 
by the graphical method we find that they correspond as 
closely as can be distinguished. The conclusion, there- 
fore, in this particular instance, is that the method 
employed has been successful. 

It should be noted that this method must be 


restricted to positive values of the objective function (Zoe 





Hyperbolic (Logarithmic) Maximization Problem 


By the 


formed 


Assume the following abstract systems 


Objective (Maximize): z 
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of logarithms the above system may be 


into the following linear system: 
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(Maximize): Z 
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By making the following substitution: 


in y = X 
ie xo = ok 


the linear system, as plotted on logarithmic coordinate 
Graphs paper adm tigurer 7 ~27.1S as /fOlLiows: 


OO Ie 2 


Ob jective (Maximize): X 5 
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IA 


Noe tn) OOK 230259 Sear ee ces os seroma (nie Sar) 
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IA 


oe ele x La (ORG foes. cme o See @ « (ie eee 
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1 P eo 9S e @ eo @ eS s » e 
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> 0.602X. de SSO2o: 3 Gabe .eeeeeiees SNeE, cet ey roa) 


HA 


0.916293 e e e e eS e 2 e e . ~(7.7.3e) 
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The inequations may be converted to equations suitable for 
submission to the simplex algorithm by the introduction of 


Slack variables, as follows: 


xX + 0.699X. = Z (Maximize) Objective) 
xy + 1.760X. + X3 =—EZaa0Zog 
xy + 1,113x, + Xp =O oS 
Xy + xX, + Xe = 1.94591 
xy + 0.602x. - Xe —lleesisiow ZS: 
xX, - Xo = TO o2ez? 
Xy - 1.474x. + Xo = 2.0/944 
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The above data, in the proper computer program format, is 
Shown in Appendix II together with the computer solution 
to the problem. The computer output indicates the solu- 
tion to this problem iss 


X, = 0.90984998- X, = 0,79135392; FUNCTIONAL = 1.46298210. 


ies my aid X5 = Ine, we take -the anti—-log=— 


rithms and find that y = 2.489° and x = 2.209” and the 


Since X 


PuneCtlOnal, = G25Lo¢ , 

In this example, as in the previous one, the results 
obtained by the graphical method (figure 7.3) compare 
favorably, insofar as distinguishable, with the results 
obtained from the computer output. The conclusions to 
be drawn, therefore, would seem to be that, in this 
particular instance anyway, the method employed has been 
successful. 

It should be noted that this method, as in the 
previous one, must be restricted to positive values of the 


objective function (2). 
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CHAPTER, Virk 
SUMMARY AND CONCLUSIONS 


An attempt has been made in this thesis to demon- 
strate the basic meaning implied in linear programming, 
and some of the ways this mathematical approach to solv- 
ing resource management problems may be applied in private 
business as well as in military supply endeavors. 

An emphasis has tried to be placed on the very 
important relationship existing between a mathematical 
description (or mathematical model) and the essence of 
what it describes, beyond any literal reality. 

It has been shown, albeit with relatively simple 
examples, that linear programming most oftentimes deals 
with rather great quantities of arithmetic in solving prob- 
lems. Problems that really need to be solved by linear 
programming, generally require the use of electronic 
digital computation equipment? otherwise, by the time some 
large problems could be solved, by hand computation, the 
problems would have become obsolete anda solution, even 
if accomplished, would be meaningless. 

One should not have received the impression that 
Minear predramming 15 a panacea: 1t is nOe.” Linear pro— 
gramming, like most everything else in life, has its dis- 


advantages and faults. 





ie 2 

Perhaps the greatest handicap (if one really could 
classify it as a handicap), to a potential user of linear 
programming, is that, almost invariably, a digital computer 
would be required to solve the problems which would actually 
be worth solving. In, and of itself this is not necessarily 
bad; it isn't bad, that is, if one always has a digital 
computer handy and the funds available with which to oper- 
ate it. It is not unheard of that when some very large 
problems have been solved by linear programming, by means 
of a digital computer, the optimum solution indicated how 
savings could be made which were so small that they were 
PWASULEIGieCne LO pay Lor Ehe time and eriore thay nad, been 
expended in solving the problem, in the first place. It 
is not the intention, here, to downgrade linear program- 
ming; on the contrary, it should be understood that linear 
programming can be, and is, a very powerful and effective 
weapon against waste and inefficiency, but, like all new 
“wonder drugs", it must be used with discretion. 

Those who are looking for an automatic decision 
making device will be disappointed with linear programming. 
To the best knowledge of this author there has not yet been 
developed a "device" or “technique" which can replace 
judgement and experience, or a “device" or technique" which 
will evaluate ethical or esthetic values and choices; at 


least, linear programming has not resolved these human 


dilemmas. 





es 

What linear programming can do, and do very well, is 
to solve certain types of mathematical models; and if the 
model has been properly constructed, and if proper data 
have been used (or, the limitations of the data are 
known), then linear programming can give an answer which 
will be the best (or the optimum) one under the particular 
set of circumstances. 

Finally, if this thesis has helped to point out a 
direction in which one may go in looking for the best way 
to manage some resource, or many resources, (especially 
resources coming within the purview of the U. S. Navy 
Supply Corps) then, the predominate purpose of this 


thesis has been accomplished. 
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APPENDIX £ 


COMPUTER PROGRAM TO DETERMINE APPROXIMATE OPTIMUM LOAD 


WEIGHT WITH MAXIMUM CARGO VALUE 


PuURvOSse 


The purpose of this program is to compute the 
approximate optimum load weight (equal to or less than a 
specified maximum load weight limit) in order to obtain 
a maximum total value of cargo for the loading problem 


as explained in Chapter VII. 


Language 


Fortran II (IBM 1620 Computer). 


Symbolic Dictlonary 


* k* 
Variable S/A I/O Descripelon 
RL S T&O Maximum allowable load (or weight) 


limit of vehicle to be loaded. 


N S T Total number of items (or pack- 
ages) to be considered for loading. 


W A at Weight of a package to be con- 
sidered for loading. 
V A i Value of a package to be consid- 
ered for loading. 
x A -- An index. Computed internally as: 
X. = v./w. 
a ane 
* 
S - Single variable; A - Array of variables. 
kK* 


c= rnp s WO Output. 





rea 


WI A O Weight of item to be loaded. 
VI A O Value of item to be loaded, 
XI A O Index of item to be loaded 


except when WI = 0.0 and 
VI = 0.0 (see below). 


WSUM S O Total weight of cargo to be 
loaded. 

VSUM S O Total value of cargo to be 
loaded. 


BEOgram (Keune 

This program utilizes the data points (representing 
weights and values of packages, or items, to be loaded into 
a vehicle having a maximum cargo weight limit) to compute 
an index (X). By ordering the data in several ways and 
performing several checking procedures, a final approxi- 
mate loading schedule is computed and is given as the out- 
put along with the total weight and value of all the pack- 
ages to be included as cargo, and the maximum allowable 
weight (objective) of the vehicle. In some instances the 
loading schedule may contain weights and values of zero 
(0.0), but indicate an index number; these will be packages 
which were in a first feasible solution, based simply on 
the index criteria, and later replaced by a package through 
subsequent checking procedures. Only packages having 
weights and values greater than zero are to be considered 


in the final loading schedule. 
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Bence jow Tech oSctLi ngs 


Sense Switch 1: When placed in the "on" position 
a listing of weights (w), values (v) and indexes (X) will 
be punched in descending order of index (X) and with 
descending order of weights (w) where two or more index 
numbers are the same. 

Sense Switch 2: When placed in the "on" position 
a listing of weights (w), values (v) and indexes (X) 


will be punched in descending order of index (X). 





Poo 
FLOW DIAGRAM FOR COMPUTER PROGRAM TO DETERMINE APPROXIMATE 
OPTIMUM LOAD WEIGHT WITH MAXIMUM CARGO VALUE 
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COMPUTER PROGRAM TO DETERMINE APPROXIMATE OPTIMUM LOAD 
WEIGHT, AND LOAD VOLUME, WITH MAXIMUM CARGO VALUE 


Purpose 


The purpose of this program is to compute the approxi- 
mate optimum load weight and load value (equal to or less 
than a specified maximum load weight limit and volume 
limit) in order to obtain a maximum total value of cargo 
for the loading problem with volume considerations as 


explained in chapter VII. 


Language 


Fortran II (IBM 1620 Computer). 


Symbolic Dict lOonary. 


* kx 
Variable S/A Pe Description 
RL S T&O Maximum allowable weight limit of 
vehicle to be loaded. 
CL S T&O Maximum allowable volume limit of 
vehicle to be loaded. 
N S I Total number of items (or packages) 
to be considered for loading. 
W A x Weight of a package to be considered 
for toading-. 
V A I Value of a package to be consid- 
ered for loading. 
Cc A I Volume of a package to be consid- 
ered for loading. 
*s -~ Single variable; A - Array of variables 


** 
tT = input: O = output 
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X A -~ An Index. Computed internally ass: 
x, = v,/(w,)(c,;) 
WI A O Weight of package to be loaded. 
VI A O Value of package to be loaded. 
er A O VOlume of package to be loaded. 
> ae A O Index of package to be loaded 


except when WI = 0.0, VI = 0.0 and 
CI = 0.0 (see below). 


WSUM Ss O Total weight of cargo to be loaded. 
VSUM S O Total value of cargo to be loaded. 
CSUM S O Total volume of cargo to be loaded. 


Program Routing 


This program utilizes the data points (representing 
weights, values, and volumes of packages, or items, to be 
loaded into a vehicle having a maximum cargo weight limit 
and a maximum cargo volume limit) to compute on index (X). 
By ordering the data in several ways and performing 
several checking procedures, a final approximate loading 
schedule is computed and is given as the output along with 
the total weight, total value and total volume of all the 
packages to be included as cargo, and the maximum allow- 
able weight (objective) and maximum allowable volume 
(objective) of the vehicle. In some instances the loading 
schedule may contain weights, values, and volumes of zero 


(0.0), but indicate an index number: these will be packages 
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which were ina first feasible solution, based, simply, on 
the index criteria, and later replaced by a package through 
subsequent checking procedures. Only packages having 
weights and values greater than zero are to be considered 


in the final loading schedule. 


penseaowlten, scetrings 


Sense Switch 1: When placed in the "on" position 
a listing of weights (w), values (v), and volumes will be 
punched in descending order of index (X) and with descend- 
ing order of weights (w) where two or more index numbers 


are the same, 
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COMPUTER PROGRAM TO GENERATE RANDOM DATA 


Purpose 


The purpose of this program is to generate simulated 
data for use in testing the approximation algorithm (i.e., 
the computer program to determine approximate optimum load 
weight with maximum cargo value). The data produced by 
this program is suitable for direct submission to either 
the above mentioned program or to the appropriate direct 


enumeration program, which follows: 


Language 


Fortran II. (IBM 1620 Computer). 


Symbolic DucuLOna ny 


* kk *& 
Variable S/A I/O Description 
NR A af Table of random numbers (1040, nine 


Gigit, fixed point, numbers). 


KK S a Number of data points (or simulated 
packages) desired as output data. 


L S aL Any random number between one (1) 
and 1039, indicating where the 
program is to begin in the random 
number table. 


% 
S - Single variable; A - Array of variables. 


** 


I - Input; O —- Output, 
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NSETS S I Number of sets of data desired, 
AX S O Simulated maximum allowable weight 
dys qy alee 
AW A O Simulated package weights. 
AV A O Simulated package values 


BEOgian RoOuUEIne 


Certain nine digit numbers are chosen from the ran- 
dom number table and through various arithmetic operations 
are converted into simulated weight and value data, as well 
as one simulated maximum allowable weight load limit. The 
First random number chosen for conversion into the simulated 
data is determined by the input variable "I". Subsequent 
numbers from the random number table are chosen by a self- 
generating random device. The data may be produced as 
output in either floating or fixed point format, or both 


(see sense switch settings below). 


Sense Switch Settings 


Sense Switch 1: When placed in the "on" position 
the output will consist of one punched card representing 
maximum allowable weight load limit, one card representing 
the number of data points to follow, and finally the 
appropriate (determined by the input variable "KK") number 
of punched cards containing both simulated weights and 


values. 


=— = 
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Sense Switch 2: When placed in the “on" position 
the output will consist of twenty-four (24) punched cards 


containing random, fixed point, numbers. 


ae ae 
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COMPUTER PROGRAMS TO DETERMINE BY DIRECT ENUMERATION, 


THE OPTIMUM LOAD WEIGHT WITH MAXIMUM CARGO VALUE 
FOR SIX, TEN AND TWELVE ITEMS (OR PACKAGES) 


EurpoOse 


The purpose of these programs is to compute all 
possible weight and value combinations, for the appropri- 
ate number of data points (6, 10, or 12), and to deter- 
mine the one set of. data points having the greatest total 
value but with a total weight equal to or less than a 
specified maximum allowable weight load limit. Since 
these following three programs are identical except for 
Varying numbers o£ input data points which are required, 
only one explanation is given here with exceptions being 


noted where necessary. 


Language 


Fortran II (IBM 1620 Computer). 


Symbolic Dieta Vv 


* * * 
Variable S/A 70 Description 
WMAX S T&0O Maximum allowable load (or weight) 


limit of vehicle to be loaded. 


NEN S ne | Dummy variable. Equivalent to 
the number of items (or packages) 
to be considered for loading. 


*s —~ Single variable; A - Array of variables 


** 
1 =| Input; O = OUTPUT. 
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SUMWS S O Total possible weight of cargo to be 
loaded. 
SUMVS S O Total maximum possible value of cargo 


to be loaded, 


AW, AV 

le ae Weights and values, respectively of 
CW, CV 

DW, DV Ss i Six packages to be considered for 

EW, EV foadi ng. 

GW, GV 

need Additional (to the above six) weights 
ne I aval tivel ie | 
RW, RV and values, respec Hood 2x2 os ten packK- 
sw, SV ages to be considered for loading. 
JR OIA: S 7" Additional (to the above ten) weights 
UW, UV and values, respectively of twelve 


packages to be considered for loading. 


Program Routine 


The maximum total value of cargo, which has a total 
weight equal to or less than the maximum allowable load 
(or weight) limit is determined by making all possible 
summations of the possible combinations of package weights 
and all possible summations of the corresponding possible 
combinations of package values. The output is in a form 
which indicates the total maximum cargo loading weight, 
the maximum allowable load (or weight) limit, and the total 
maximum cargo value. In addition the output contains a 
loading schedule having a heading (of alphabetic letters 


A, B, C, etc.) corresponding to the individual data points 
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and either the number 1 or the number 2 underneath the 
respective alphabetic letter. The number 2 indicates 
that the package is to be loaded and is included in the 


loading schedule. 
Sense Switch settings 


sense Switen 2s When placed in/the~ “on” position 
Pe cOUrE DUT TCOnStSting Of Cach individual calcularicn wi! 
be produced. The output will contain the combinations of 
packages being considered and the total combined weight 
amdevalue Of Ehat. computation 

wemse SWiECn 3: “When placed Am Ene: “on” spesteso, 
the program will pause at the completion of all computa- 
tions. By pressing "start", on the computer console, the 
program will continue with another set of data. When sense 
SWwiten) 3 1S im the “OLE” position the program wills continue 
with another set of data automatically, when completed with 


the previous computations. 





MbOw IlLAGRAMe FOR COMPUIER PROGRAM "lO DETERMINE, 
BNUMB RATION. “THE OPTIMUM MOAD WEIGHT Wire 
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